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abstract An extensive spectrum of light non-strange qq states up to a mass 
of 2400 MeV has emerged from Crystal Barrel and PS172 data on pp —>■ 
Resonance —>■ A + B in 17 final states. These data are reviewed with detailed 
comments on the status of each resonance. For I = 0, C = +1, the spectrum 
is complete and very secure. 

O 

O ■ Six 'extra' states are identified. Four of them have / = 0, C = +1 and spin- 

parities predicted for glueballs. Their mass ratios agree closely with predictions 
Q . from Lattice QCD calculations. However, branching ratios for decays are not 

Q [ flavour blind and require mixing with neighbouring qq states. This mixing 

l/^ I remains controversial. 

Two further states, 772(1870) and 7r2(1880), are candidates for a pair of hybrids. 

> ; Their decay branching ratios BR[T]2il870) /2(1270)r/]/5i?[r/2(1870) 

1^ ; a2(1320)7r] and 5i?[7r2(1880) ^ a2(1320)r/]/5i?[7r2(1880) ^ /2(1270)7r] are 

■ much larger than predicted by SU(3) for qq states. 

(N 

^ ^ A fourth category of anomalously light 0+ mesons is made of (t(540), k(750), 

O ' ao(980) and /o(980), forming a non-gg nonet. Their properties are reviewed 

■ using data from both Crystal Barrel and elsewhere. 
CD ■ 
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1 Introduction and Survey of Results 

Many new resonances have appeared in Crystal Barrel and PS172 data on pp 
annihilation in flight. They follow a remarkably simple pattern. The objective 
of this review is to outline techniques used in the analysis and also to assign a 
'star rating' to each state. Hopefully the star rating will assist comparison with 
model predictons, where it is important to know which states are identified 
most reliably. Techniques of analysis may be of interest to further experiments 
which will produce overlapping information. Results from 35 journal publica- 
tions in collaboration with the St. Petersburg group are summarised. 

The main idea is to study s-channel resonances: 

pp — > Resonance A + B, 

by full partial wave analysis of 17 final states A + B listed in Table 1. Cou- 
pled channel analyses have been carried out, following the classic techniques 
applied to the nN system in 1960-70, imposing the fundamental constraints 
of analyticity and unitarity. 
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Table 1 

Channels studied; decays 77 — 77 and r] — 37r° are both used, and 77' — •qir'^'K^; also 
ui — > 7r°7 and u — ^ tt'^ it" ir^ . 

The majority of final states are made up only of neutral particles. The Crystal 
Barrel has excellent 7 detection in a barrel of 1380 Csl crystals with good 
energy and angular resolution (see Section 2). Crystals cover 98% of the solid 
angle. Detection of charged particles is limited to 71% of the solid angle. 
Neutral final states fall cleanly into four categories of isospin I and change 
conjugation C, as shown in Table 1; 7r~7r+ contains both 7 = and 1. 

The largest number of decay channels is available for 7 = 0, C = +1. For the 
Ti~Ti^ final state, da/dQ and polarisation data are included from earlier ex- 
periments [1,2]. The polarisation data play a vital role in identifying helicities 
and hence separating pp states from ^P2- There is also a second point: 

oc Im{FlF,) (1) 

where Fi and F2 are helicity combinations discussed later. Differential cross 
sections depend on moduli squared of amplitudes and real parts of interfer- 
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ences. The combination of differential cross sections and polarisations deter- 
mines both magnitudes and phases of amplitudes uniquely. In any one data 
set, significance levels of resonances are similar to those of vrA^ analyses above 
1700 MeV; for / = 0, C = +1, the analysis of 8 channels with consistent 
parameters leads to extremely secure amplitudes and a complete set of the 
expected qq states in the mass range 1910-2410 MeV, each with a statistical 
significance > 25(7. 

The spectrum shown in Fig. 1 [3] is remarkably simple. Almost within errors, 
observed states he on parallel trajectories of (where M is mass) v. radial 
excitation number n. The mean slope is 1.14 ± 0.013 GeV^ per unit of n. 

These trajectories resemble the well-known Regge trajectories; however on 
Regge trajectories, states alternate in / and C (for example pi, /2, ps, /4, p^). 
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Fig. 1. Trajectories of / = 0, C = +1 states. Numbers indicate masses in MeV. 

Similar spectra are observed for other / and C, but are less complete and 
less precise. Fig. 2 shows the spectrum for / = 1, C = +1 [4,5]. Since quarks 
arc nearly massless, one expects results to be almost degenerate with / = 0, 
C = +1. This is indeed the case. However, the separation between and 
is much poorer for 7 = 1, C = +1 because no polarisation data are available. 

For J = 1, C = — 1, the observed spectrum is shown in Fig. 3(b) [6,7]. Here 
polarisation data for pp — > tt~7t~^ are available and do help enormously in 
separating magnitudes and phases of amplitudes. As a result, the observed 
spectrum is complete except for one high-lying J^'^ = 1 state and the ^Gs 
pp state. The ^D^ states are particularly well identified. 
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Fig. 2. Trajectories of / = 1, C = +1 states. 




Fig. 3. Trajectories of C = —1 states with (a) and (b) / = 0, (c) and (d) I = 1. In 
(b), the ^D^ trajectory is moved one place right for clarity; in (d), is moved 
one place left. 

Figs. 3(a) and (b) show the spectrum of / = 0, C = — 1 states [7]. Errors are 
the largest here. Statistics for the Lorj channel are typically a factor 6 weaker 
than for a;7r°. The second available channel is a;7r°7r° and is difficult to analyse 
because of broad 0"*" structure in 7r°; those states which are observed are again 
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almost degenerate with I — 1, C — —1. 

All observed resonances cluster into fairly narrow mass ranges (i) 1590-1700 
MeV, (ii) 1930-2100 MeV, (iii) 2240-2340 MeV. These 'towers' of resonances 
are illustrated in Fig. 4. 
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Fig. 4. Clustering of resonances in mass. 

Sections 1-9 review qq states in detail. There are four 'extra' states with 
/ = 0, C = +1 having masses close to Lattice QCD predictions for glucballs 
and the same J^. Section 10 reviews these glueball candidates and the mixing 
of /o(1500) with the nearby nonet made of /o(1370), ao(1450), i^o(1430) and 
/o(1710). Section 11 reviews the light scalar mesons: cr(540), fi;(750), /o(980) 
and ao(980). Section 12 take a brief look at what is needed in future. 



2 The Crystal Barrel experiment 

This experiment ran at LEAR, the p storage ring at CERN, using nine p beam 
momenta of 600, 900, 1050, 1200, 1350, 1525, 1642, 1800 and 1940 MeV/c. The 
corresponding pp mass range is from 1962 to 2409 MeV. Most of the data were 
taken during the last 4 months of LEAR operation, August - December 1996. 
Data from the PS172 experiment wth a polarised target will also be used. 
This experiment ran during 1986 from 360 to 1550 MeV/c; the extension of 
the mass range down to 1910 MeV is important. 

The Crystal Barrel detector is illustrated in Fig. 5. Details are to be found 
in a full technical report [8]. A 4.4 cm liquid hydrogen target was surrounded 
sequentially by multiwire chambers for triggering (labelled 5 on the figure), a 
jet drift chamber (4) for detection of charged particles, a barrel (3) of 1380 
Csl crystals covering 98% of An solid angle, and a magnet (1) producing a 
solenoidal field of 1.5T (though this was run at l.OT for all-neutral data). 
For running after 1994, the multiwire chambers were replaced by a silicon 
vertex detector. However, for detection of neutral final states, the chambers 
and silicon vertex detector acted simply as vetos over a sohd angle 98% of An. 
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Fig. 5. Overall layout of the Crystal Barrel detector, showing (1) magnet yoke, 
(2) magnet coils, (3) Csl barrel, (4) jet drift chamber, (5) multi-wire proportional 
chambers, (6) liquid hydrogen target. 

The Csl crystals measured photon showers with an angular resolution of 20 
mrad in both polar and azimuthal angles. The photon detection efficiency was 
high down to energies below 20 MeV. The energy resolution was given by 
/\E/E = 0.025/£'^/^, with E in GeV. In order to reject events which failed 
to conserve energy, the total energy deposited in the Csl crystals was derived 
on-line; the trigger selected events having a total energy within 200 MeV of 
that in pp interactions [9] . 

Incident p were detected by a fast coincidence between a small proportional 
counter P and a Si counter Sic of radius 5 mm at the end of the beam-pipe and 
only a few cm from the entrance of the liquid hydrogen target. Four further 
silicon detectors located at the same point and segmented into quadrants were 
used to steer and focus the beam. Two veto counters were positioned 20 cm 
downstream of the target and inside the Csl barrel; they provided a trigger for 
interactions in the target. The data-taking rate was ~ 60 events/s at a beam 
intensity of 2 x 10^ p/s. 

One gets a general impression of the cleanliness with which neutral particles 
were detected from Fig. 6. Panels (a) and (b) show rj peaks in the t]-)] final 
state, from 27 events in (a) and from 37r° events in (b) [10]. As an illustration 
of selection procedures, consider pp rjr]. Cuts were applied on confidence 
levels (CL) as follows: (i) CL{7r^-f-f) to remove vrOvrO, (ii) (^^(vr^) < 0.01%, 
(iii) CL{riri) > CL{riri'), (iv) CL{riri) > 0.1, In the few cases where there 
was more than one 7777 solution, selected events were required to have a CL 
at least a factor 10 better than the second. The entire process was followed 
through a Monte Carlo simulation discussed below. Observed backgrounds 
agreed closely with this simulation after allowing for branching ratios of the 
background channels. The t] peaks in both (a) and (b) were centred (purely 
from kinematics of the production reaction) within 1-2 MeV of the Particle 
Data Group value (PDG) [11] at all beam momenta; this checks systematics of 
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Fig. 6. Distributions of (a) M(77) near the r] from events fitting 7777, after cuts; 
(b) M{3tt^) near the 77 from ■qSir'^ events; (c) M(77) near the rf' from 7777 data; (d) 
M{r)Tr^Tr'^) near the 77' from r/r/Tr^Tr^ events. Data are at 1800 MeV/c. 

the energy calibration. Fig. 7 shows angular distributions derived from 47 and 
87 events, independently normalised to Monte Carlo simulations. The close 
agreement in both shapes and absolute magnitudes is a major cross-check. 

Figs. 6(c) and (d) show rj' peaks from 2j and rjir^TT^ decays. Statistics are 
lower than for 7777, but angular distributions and normalisations agree between 
47 and 87 for pp — > 77'77 within statistics [10], and also for pp — > 77^° [12]. 

Showers were selected from adjoining groups of nine Csl crystals. The relation 
between shower energy and pulse-height in crystals was obtained [8] by fitting 
the 7r° mass for 77 events, using all crystal combinations. Three variants of this 
calibration were tried. The primary effects of different calibrations are removed 
by the overall kinematic fit to energy-momentum conservation; this uses the 
beam momentum, known to better than 0.1%. The result is that different 
cahbrations altered < 1% of selected events. However, the performance was 
better than that. Between one calibration and another, the overall change 
in accepted events was a factor 10 smaller; i.e. changing from one selection 
to another, the few events which were lost were compensated by an almost 
equal number entering the sample. A more important consideration was the 
identification of crystals which malfunctioned. For most running conditions, 
there were at most 1 or 2 crystals which failed in some way (noisy or dead). 
These were identified and cut completely from both the data and the Monte 
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Fig. 7. Comparison of angular distributions for r/r/ from 47 events (black circles) and 
87 (open triangles). Full curves show the resTilt of the amplitude analysis. Dotted 
curves illustrate the acceptance for three momenta; small undulations are due crystal 
sizes and overlap of showers between crystals, well simulated by the Monte Carlo. 

Carlo simulation. Otherwise a faulty crystal could lead to a dip in the angular 
distribution over the angles covered by 9 adjoining crystals. 



2.1 Rate dependence 



We now consider the absolute normalisation of cross sections. This is impor- 
tant in an amplitude analysis covering many momenta. A potentially serious 
rate dependence was found, but ultimately affected normalisations by only 
~ 3 — 6%, varying slowly with beam momentum. It arose from pile- up in the 
Csl crystals. Beam intensities were 2 x 10^/s and interaction rates in the target 
were typically 6 x 10^/s. 

The scintillation from Csl has two components, one with a risetime of 6/is, 
but with a fall-time of 100/xs. The discriminators used pole- zero cancellation 
which was supposed to eliminate the long component. But in practice large 
energy deposits in one crystal could produce some ringing in its discriminator 
output. This led to significant pile-up between one event and the next. A 
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Fig. 8. Rate dependence observed at (a) 1800, (b) 1525 (c) 1350 and (d) 1050 MeV/c 
for the sum of (4 + 5 + 6 + 7 + 8)7 events. Lines show the fitted rate dependence. 

detailed discussion is to be found in Ref. [10]. 

The normalisation was derived from beam counts P. Sic, the length and den- 
sity of the target, the observed number of events and the Monte Carlo simula- 
tion of detection efficiency. The general off-line procedure was to reconstruct 
45 exclusive final states, each one kinematically fitted to channels with 4 to 
IO7. Rather few have 9 or 10 7. Fig. 8 shows as a function of beam rate the 
total of events fitting with confidence levels > 10% to any channel with 4 
to 8 7. Typically ~ 20% of triggers fit such channels. There is obviously a 
strong rate dependence and it was necessary to extrapolate cross sections to 
zero beam rate. The rate dependence varied with beam steering. At several 
momenta, separate running conditions gave different slopes, but extrapolated 
values were consistent within 3%. 

Fortunately, the rate-dependence is found to kill genuine events of all types in 
the same way. Let N{n'y) be the number of events fitting above 10% confidence 
level with n photons, and let A^(4 — 8)7 be the sum of these from n = 4 to 
8. Fig. 9 displays ratios N{n'~f)/N{4 — 8)7 for 47, 67, 87 and for tt^tt^ as 
a function of beam rate. Any variation with rate is < 3%. Neither is there 
any observed change with beam rate in 7r°7r° angular distributions or Dalitz 
plots of 67 events [13]. The conclusion is that angular distributions can be 
extracted safely at any beam rate; but the normalisation must be obtained by 



10 





Fig. 9. Fractions of (a) 47, (b) 67, (c) 87 events and (d) vr'^vr*^ events, out of the 
total of final (4 + 5 + 6 + 7 + 8)7 events at 1800 MeV/c v. incident p rate. 

extrapolating to zero beam-rate. A cross-check is also satisfied between n^n^ 
data from Crystal Barrel and symmetrised 7:~tt~^ differential cross sections 
from PS 172, see below. 



2.2 Monte Carlo and Data Selection 



Data arc fitted kincmatically to a large number of physics channels: 45 for 
(4 — 10)7. In order to assess branching ratios to every channel and cross-talk 
between them, at least 20,000 Monte Carlo events are generated for every one 
of the 45 fitted channels, using GEANT. In the first approximation, events 
fitting the correct channel determine the reconstruction efficiency in that 
channel. Events fitting the wrong channel measure the probability of cross- 
talk Xij between channels i and j. More exactly, a set of 45 x 45 simultaneous 
equations is solved containing on the left-hand side the observed number of 
fitted data events Di, and on the right-hand side reconstruction efficiencies 
and true numbers of created events Ni in every channel plus terms allowing 
for cross-talk between channels: 

A = eiNi + J2 XijNj. (2) 

The solution is constrained so that the numbers of real events in every 
channel are positive or zero. 
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This procedure is carried out for several confidence levels (1, 5, 10, 20%) and 
using a wide variety of selection procedures. A choice is then made, optimising 
the ratio of signal to background. It turns out that this ratio is not very 
sensitive to confidence level over the range 5-20%. Hence selection criteria are 
not usually critical. 



3 How to establish resonances 

One objective of this review is to assign a star rating to all resonances. This 
is to some extent subjective, but some assessment is better than none, as a 
guide to further work. Before coming to details, it is necessary to review the 
general approach to identifying resonances. 

It is sometimes argued that magnitudes and phases of all amplitudes should 

be fitted freely, so as to dcmonstrac that phases follow the variation expected 
for resonances. Though that is a fine ideal, it is unrealistic and leads to prob- 
lems. Phases can only be determined with respect to some known amplitude. 
As examples, phase variation of the Zq and the Jj^ may be found from in- 
terferences with electrweak amplitudes. In i:^]) scattering at low energies, the 
A(1232) is established as a resonance by the fact that its cross section reaches 
the unitarity limit; interferences with S-waves then establish its resonant phase 
variation, since there is not enough cross section for S-waves to resonate too. 
Further examples are rather rare. 

The usual assumption is that isolated, well defined peaks far from thresholds 
are assumed to be resonances. The 7r2(1670) and /4(2050) are examples. The 
phases of other amplitudes may be determined using these and other estab- 
lished resonances as interferometers. This is the procedure followed in Crystal 
Barrel analyses. 

The /4(2050) appears strongly in data for / = 0, C = -|-1. It acts as an 
interferometer and determines phases of other partial waves. They have a 
similar phase variation with s and must therefore also resonate. In pp — > tt^tt"*", 
differential cross sections establish real parts of interferences and polarisation 
data measure the imaginary parts of the same interferences. Together, these 
data establish quite precisely the phase variation of all partial waves with 
respect to /4(2050). Further strong peaks appear as /4(2300) and 772(226?). 
They act as further interferometers, loosely coupled themselves in phase to 
the tail of /4(2050). 

A strong p3(1982) appears in pp — > 7r~7r"^. In turn, this acts as an interferom- 
eter for 7 = 1, C = — 1. The amplitude analysis grows like a rolling snowball. 
The phases of other partial waves at low mass are measured with respect to 



12 



P3(1982) and establish resonances in other partial waves at low mass. Vari- 
ations of magnitudes are consistent with the observed phases. Peaks appear 
also in the mass range 2200-2300 MeV, all correlated in phase. Collectively, 
their phases are loosely coupled to the tails of lower mass states. 

For I = 1, C = +1 and J = 0, C = — 1 there are no polarisation data. 
Phase information from differential cross sections arise only via real parts of 
interferences. When relative phases are close to or 180°, errors of phases are 
large. Without some smoothing of the phases to follow resonant behaviour (or 
some other analytic form) , noise in the analysis can run out of control and can 
lead to a forest of narrow peaks. Such a result is unphysical, since amplitudes 
should be analytic functions of s. Smoothing to analytic forms improves the 
determination of partial wave amplitudes dramatically. 

The technique used for ttN elastic scattering is to constrain the analysis using 
fixed t dispersion relations. This works successfully up to 800 MeV. Above 
that energy, amplitudes are parametrised directly with analytic forms. 

The approach used in Crystal Barrel analysis is to use analyticity as a foun- 
dation, parametrising amplitudes from the outset in terms of simple Breit- 
Wigner resonances of constant width, or constants, or linear functions of s. 
Even then, caution is necessary. If a linear function develops a phase variation 
> 90°, it is likely to represent a resonance. It is well known to be dangerous to 
fit data with amplitudes having similar characteristics; orthogonal functions 
are needed. In reality, one can usually replace linear functions having rising 
phases as the tails of resonances well below the pp threshold, e.g. 7r(1800) or 
P3(1690); their phase variations above the pp threshold are small and under 
reasonable control. Linear variations with falling phases are a possibility but 
in practice are not observed. A valuable constraint in the analysis is that am- 
plitudes for all partial waves except ^5*0 and ^5*1 go to zero at the pp threshold; 
these two go to scattering lengths and are the most difficulty to deal with. 

There is a practical problem with the pp threshold. Fig. 10 shows integrated 
cross sections for two-body channels. For neutral final states 7r°7r° and rjrj, 
only = 0+, 2"*", 4"*" and G"*" are allowed by Bose statistics. There are no 
contributions from pp S-waves, which have quantum numbers J^*" = O""*" {pp 
^Sq) and 1 {pp ^Si). One sees in Fig. 10 that neutral cross sections go to low 
values near the pp threshold. However, the 7r~7r+ cross section goes to large 
values, following the familar 1/v law. The pp total cross sections follows a 
similar behaviour. The S-wave amplitudes near threshold need to contain this 
feature, and it becomes difficult to separate the 1/v behaviour from resonances 
at the pp threshold. 

A further problem is that singlet pp states do not interfere with /4(2050) (the 
triplet state in pp). 
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Fig. 10. Integrated cross sections for 2-body scattering, compared with the partial 
wave fit. In (a), black circles are data of Hasan et al. and open squares are those 
of Eisenhandler et al. In (b), (c) and (d) the integration is from cos 9 = to 0.875 
and is over only one hemisphere; the r]r]' data are averaged over 47 and 87 data. 
For 7r~7r"'", the integration is over the full angular range. In (c), the dashed curve 
shows the O"*" intensity. 

An experimental problem is that Crystal Barrel data at the lowest momenta 
600 and 900 MeV/c straddle the lowest tower of resonances. These momenta 
correspond to masses of 1962 and 2049 MeV. In several cases, what can be said 
is that there are definite peaks, inconsistent with the tails of resonances be- 
low threshold. However, because of the limited number of available momenta, 
masses and widths become strongly correlated; precise masses are subject to 
some systematic uncertainty, which is correlated over resonances of different 
via interference effects. Each of these cases needs to be examined critically. 
For 7 = 0, C = +1 and 7 = 1, C = -1, the PS172 data for pp n-n+ down 
to 360 MeV/c are very important in extending the mass range to 1912 MeV 
in small steps of p momentum. 

3. 1 Star Ratings of Resonances 

Resonances will be given ratings up to 4*. The highest class requires observa- 
tion of 3 or more strong, unmistakable peaks and a good mass determination 
(with error SM < 40 MeV). Such states arc equivalent to those in the sum- 
mary table of the Particle Data Group - better than some cases like ao(1450). 
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which until 2003 was observed only in a single channel. 

States are demoted to lower star ratings if the mass determination is poor, 
particularly if F > 300 MeV. The 3* resonances are reasonably well estab- 
lished, usually in two strong channels with errors for masses < ±40 MeV; 
some are established in 3 channels but with 6M in the range 40-70 MeV. 

The 2* resonances need confirmation elsewhere. They arc observed either (i) 
in one strong channel with 5M < 40 MeV or (ii) in 2 strong channels with 
sizable error in the mass. 

Finally, 1* states are tentative: weak channels and poor mass determination. 



4 The Partial Wave Analysis 

The foundation for the analysis was laid in 1992-4 while fitting PS 172 data on 
pp Tr^TT"*" [1]. Early in this analysis it became clear that there is a large J — A 
amplitude which can be identified with the well known /4(2050). At momenta 
> 1 GeV/c, a prominent feature of the polarisation data, illustrated below 
in Fig. 13, is a polarisation close to +1 over much of the angular range and 
over most beam momenta, broken by some rapid fluctuations. If there were a 
single /4 resonance and non-resonant amplitudes in other partial waves, the 
polarisation would necessarily show a dramatic energy variation through the 
resonance. This is in complete contradiction with the data. The only way to 
fit the data successfully is to have resonances in many partial waves. 

If one invokes analyticity to relate real and imaginary parts of amplitudes, the 
two available sets of data, (polarisations and differential cross sections), can 
in principle lead to a unique set of partial waves. The first steps in analysing 
Crystal Barrel data were to add high statistics data for pp — > tt^tt", rjrj, rjrj' and 
rjTT^TT^, all having 7 = 0,C = -|-1. It rapidly became apparent that all these 
data could likewise be fitted in terms of resonances and simple backgrounds. A 
coupled channel analysis led to consistency between all sets of data and rather 
precise determinations of resonance parameters [13-15]. The final analysis for 
these quantum numbers is in Ref. [3]. 
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4-1 Formulae for Amplitudes and Resonances 

It is easiest to condider one reaction as an example: pp — > ao(980)7r° — > r;7r°7r°. 
The amplitude for one partial wave is written 

Ppp ge^^B,{p)BL{q) . . , 

^ - -J M^-s-iMV ^^y'^^^^'^ + F[s,,)Z,,\. (3) 

Here gi is a real coupling constant, and the strong- interaction phase for this 
partial wave. Blatt-Weisskopf centrifugal barrier factors B are included explic- 
itly for production with orbital angular momentum £ from pp and for coupling 
with angular momentum L to the decay channel (momentum q). The first 
centrifugal barrier includes the correct threshold dependence on momentum p 
in the pp channel. Formulae for B are given by Bugg, Sarantsev and Zou [16]. 
The radius of the centrifugal barrier is optimised at 0.83 fm in Ref . [3] . It is 
important only for the high partial waves where the barrier is strong and well 
determined. For L = 4 waves, the barrier radius is increased to 1.1 fm. The 
factor 1/p in eqn. (3) allows for the flux in the pp entrance channel, and ppp 
is the phase space for this channel, 2p/ ^/s. Near the pp threshold, the S-wave 
cross section is then proportional to 1/p^/s, the well-known 1/v law. 

Strictly speaking, the Breit-Wigner amplitude should take the general form 

/ = l/[M2-s-m(s) -iMr(s)], (4) 

where r(s) is the sum of widths to decay channels. For each channel, T{s) is 
proportional to the available phase space. This is the well known Flatte form 
[17], which will be used below for special cases such as /o(980) and a0(980); 
these resonances lie close to the KK threshold and the s-dependence of the 
KK width is important. There is a dispersive term m{s) given by Tornqvist 
[18]. For present purposes, these refinements arc not relevant. The analysis 
concerns high masses where many decay channels arc open. Many of the decay 
channels are unknown. The approximation is therefore that F can be taken 
as a constant; m{s) is then zero. Backgrounds are usually parametrised as 
constants or as resonances below the pp threshold. 

For the 777r°7r° final state, there are two possible amplitudes for production 
of ao(980) in r^vr; F{syi) is the Flatte formula for the ao resonance between 
particles 1 and 2, and F{sri) is the corresponding amplitude between particles 
1 and 3. The functions Z12 and Z13 describe the dependence of the production 
amplitude on the production angle 6 in the centre of mass and on decay angles 
a, (5 for decay in the resonance rest frame. These functions are written in terms 
of Lorentz invariant tensors, following the methods described by Chung [19]. 

Partial waves with = 2+ and 4+ couple to pp with i = J ± 1. Multiple 
scattering through the resonance is expected to lead to approximately the 
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same phase for both £ values. The ratio of coupling constants gj+i/gj^i is 
therefore fitted to a real ratio rj; this assumption can be checked in several 
cases and data are consistent with this assumption. 

A general comment concerns possible i-channel contributions. Singularities 
due to those exchanges are distant and therefore contribute mostly to low 
partial waves. The high partial waves with J > 2 are consistent with no 
background. Secondly, t-channel exchange amplitudes acquire the phases of 
resonances to which they contribute and may be considered as driving forces 
for those resonances; an illustration is the well-known Chew-Low theory [20], 
where nucleon exchange affects the line-shape of A(1232). In the present anal- 
ysis, line-shapes of resonances cannot be distinguished in most cases from 
simple Brief- Wigner resonances with constant widths. 

4-2 Illustration of spectroscopic notation 

As an example of spectroscopic notation, consider pp — > cutt. The initial pp 
state may have spin s = or 1. For singlet (s = 0) states, the total angular 
momentum is J = i. For triplet states {s = 1), J = i or i ± 1. The parity is 
P = (—1)^ and C = (—1)^"'"''. In the uutt exit channel, the spin s is 1 for the uu. 
As examples, J^*^ = 1 may couple to initial pp ^Si and ^Di partial waves 
having ^ = or 2; the exit channel has L = 1. In the analysis, ^5"! and ^Di 
partial waves are fitted with a ratio of coupling constants rj = g^^j^i/ gi^=j_i. 
A second example is J^*" = 3^ . This couples to pp singlet states and decays 
to LJTT^ with L = 2 or 4. These two partial waves arc again fitted with a ratio 
of coupling constants rj = gL=j+i/gL=j-i- Multiple scattering through the 
resonance is expected to lead to approximately the same phase for both 
values of i or L. It is indeed found that all phase differences between £ = J ± 1 
lie within the range to ±15°. Where possible, ratios rj are taken to be real, 
in order to stiffen the analysis. 

The ujriTi^ data arc fitted to sequential 2-body processes: a;a2(1320), cjao(980), 
6i(1235)?7 and C(j(1650)7r*^. As regards notation, consider uja2. Spins 1 of the u 
and 2 for the a2 combine to total spin s' = 1, 2 or 3. The initial pp ^F^ states 
with J^'-^ = 3"*"" may decay to a; 02 with L = 1, 3 or 5. In practice, the lowest 
L value is always dominant; in every case L values above the second may be 
omitted because of the strong centrifugal barrier. It is however necessary to 
consider all possible values of total spin s'. 
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4-3 General comments on procedures 



The selection of amplitudes is kept as simple as possible, consistent with the 
data; introducing amplitudes whose significance is marginal introduces noise. 
A rule-of-thumb is that one can normally remove amplitudes which improve 
log likelihood by < 20. For / = 0, C = +1 where 8 data sets are available 
and the solution is strongly over-constrained, amplitudes contributing < 40 
are omitted; they contribute < 0.3% of the total intensity in all cases. 

Programmes can optimise M and F as well as coupling constants and the 
barrier radius R. In practice it is important to scan all M and F through at 
least 7 steps (freeing other M and F) so as to check that there is a well defined 
parabolic optimum in log likelihood. Movements of other M and F reveal 
correlations. For a few poorly determined resonances near the pp threshold, log 
likelihood flattens out for masses below threshold. Likewise, poorly determined 
widths may tend to flatten off or run away at large values; this can arise via 
cross-talk with other resonances having the same quantum numbers at higher 
or lower masses. Resonances with poor behaviour of this type are assigned a 
low star rating or are absorbed into background terms where possible. 

It is bad practice to introduce broad isotropic backgrounds. These interfere 
over the whole Dalitz plot and can perturb resonances. A symptom of instabil- 
ity is the development of large destructive interferences between components. 
This is checked by printing out the integrated contributions from every am- 
plitude and every pair of interferences. 

Well known resonances such as /2(1270), 02(1320), ao(980), 6i(1235), etc are 
flxed at masses and widths quoted by the PDG; the D/S ratio of amplitudes 
for decay of 6i(1235) (and other resonances) are fixed at PDG values. 

A general point is that programmes must be as fast and as simple as possible, 
since one needs to run hundreds of variants of the fit with varying sets of 
amplitudes. A simple though laborious procedure to look for alternative solu- 
tions is to flip the signs of the smaller amplitudes one by one and in pairs and 
check whether the solution returns to its starting point or a different solution. 
Signs of r parameters are also flipped. In practice Barrelet ambiguities do not 
survive in the coupled-channel fits. Only one ambiguous solution is found, for 
rjTi^ as described below. The bigger problem arises with poor solutions where 
several amplitudes are strongly correlated. The addition of each new set of 
data improves the stability of the solution greatly; one new channel is worth 
a factor 4 statistics in existing data. Generally speaking, 50,000 events per 
channel are sufficient; further statistics do not improve the solution and one 
has to beware that errors become systematic rather than statistical. 

Cross-talk between high partial waves arises because of hmitations on accep- 
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tance for small production angles. For Crystal Barrel data, cross-talk begins 
to develop between i = 4 and £ = 5 partial waves. A general guide is that it is 
necessary to have an acceptance which extends to halfway between cos^^ = ±1 
and the first zero of the highest Legendre polynomial. To identify i = 5 partial 
waves requires acceptance to cos 9 — ±0.96 or ^ < 6°, where 9 is the centre of 
mass production angle. In practice this is available in the backward hemisphere 
because of the Lorentz boost, but is not available in the forward hemisphere. 
An obvious point is that bins where acceptance varies rapidly with scatter- 
ing angle should be removed or flagged. For all-neutral flnal states, angular 
distributions are symmetric forward-backward because of conservation of C- 
parity; physically, neutral products do not distinguish between initial p and 
p. A feature of the PS172 data is that the angular acceptance is complete to 
cos^ — ±1 since the beam went through the multi-wire detectors. 

Adjoining resonances which overlap by more than their half-width are poorly 
separated. This makes separation of and ^Di resonances difficult, since 
they alternate at steps of ~ 150 MeV and have widths typically 250 MeV. 

One sometimes finds that resonance parameters are determined sensitively in 
simple channels, e.g. ujt), despite low statistics. Higher statistics may exist in 
more complex channels, e.g. ujir^n^. However, analysis of these more complex 
channels may give poorer sensitivity, for example due to the broad a in the 
ua channel. [Here a denotes the whole tttt S-wave amplitude]. In this case it 
may be desirable to fix resonance parameters for these more complex channels. 
There is an alternative technique for dealing with this problem. It may also be 
used when searching for weak decay modes of resonances. A penalty function 
is introduced. This penalty function takes the form of contributions to x^'- 

5A \i' 

Here A is a parameter and Aq a target value. The denominator 5A is adjusted 
so that each questionable ingredient contributes = 9 (i.e. 3a). Those am- 
plitudes which are really needed feel little influence from this small penalty 
function compared to the influence of data points; amplitudes which are unsta- 
ble or not needed are well constrained. This is particularly useful for searching 
for small branching ratios; here Aq is set to zero. This procedure is applied 
for I — 1, C — +1, where some instabilities are discussed below; the penalty 
function contributes ~ 650 to compared with ~ 650, 000 for log likelihood. 
Incidentally, the program optimises both ioT some data sets (e.g. binned 
2-body angular distributions) and log likelihood (for 3-body channels). 

There is a useful trick in checking of resonances. It is difficult to exhibit 
the angular dependence of some amplitudes, for example high a27r; the 

angular dependence Z of eqn. (3) involves complicated correlations between 
production angle 9 and decay angles a, (3 of the a2. Putting any new resonance 
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into the analysis always gives some improvement. Some of this improvement 
comes from the correct line-shape and some from the correct angular depen- 
dence. It is useful to substitute different (wrong) angles and examine how 
sensitive the analysis is to the right choice of angles and correlations. 



5 Review of 7 = 0, C — +1 resonances 



Most resonance from the Crystal Barrel in-fiight analysis are listed by the 
Particle Data Group in a section called 'Other Light Unfiavoured Mesons', pp 
010001-539 of Ref. [11]. The reason is that the Particle Data Group requires 
confirmation from a separate experiment, rather than from different channels 
of data in one experiment. There is very little other coverage of this mass range 
from other experiments. Some results are combined by the PDG with those 
of other groups in the usual way, some are not. Many of the resonances de- 
scribed here are observed in several decay channels and are better established 
that other entries from single groups, so this creates a confusing situation. 
Resonance parameters from the final / = 0, C = -|-1 coupled-channel analysis 
[3] are given in Table 2. 



5.1 2-body channels 



Resonances with — 0~^, 2"*" and 4"*" arc observed in 2-body channels, where 
polarisation data are available. The first step was to determine these ampli- 
tudes. The extensive differential cross section and polarisation data come from 
Eisenhandler et al. [2] and PS172, Hasan et al. [1]. They are compared with 
the fit in Figs. 11-13. 

Why are polarisation data needed? For pp n^n^, there are only triplet 
amplitudes with J = L± 1: ^^i, ^Pq, ^Di, ^D^, .... It is convenient to use two 
amplitudes which measure transfer of helicity -1-1/2 of the proton in the initial 
state to ±1/2 in the final state. These are denoted F++ (hehcity non-flip) and 

F+_ (helicity flip). [F = and F_+ = In terms of partial wave 

amplitudes and Legendre polynomials, 

F^^ = Wi^J + l)fi^Pj{cose) (5) 
^ J 

= i E -^^^fi+PKcosO), (6) 



J 



V2jTlf++ = y/jTj_i,j - y/jTlTj+^,j (7) 
v/2JTT/+- = vOTTrj-i,j + ^Tj+,,j. (8) 
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J^^ * Mass M Width T r A5 Observed 







(MeV) 


(MeV) 




(ryvrvr) 


channels 


6++ 


- 


2485 ± 40 


410 ± 90 





245 




4++ 


4* 


2283 ± 17 


310 ± 25 


2.7 lb 0.5 


2185 


7r7r,7?r?,/2T/,VES 


4++ 


4* 


2018 ± 6 


182 lb 7 


0.0 lb 0.04 


1607 


7r7r,7777, r/r/', a27r,(/2r7),PDG 


4-+ 


2* 


2328 ± 38 


240 lb 90 




558 


ao7r,(a27r) 


3++ 


4* 


2303 ± 15 


214 lb 29 




1173 


[/2r/]L=i,3,/2r/' 


3++ 


4* 


2048 ± 8 


213 lb 34 




1345 


[/2'?]L=i,3,a27r,ao7r 


2++ 


4* 


2293 ± 13 


216 lb 37 


-2.2 lb 0.6 


1557 


TnT,r]r],rir]'j2r],a27r 


2++ 


4* 


2240 ± 15 


241 ± 30 


0.46 lb 0.09 


468 




2++ 


4* 


2001 ± 10 


312 ± 32 


5.0 lb 0.5 


168 


TTTT, r]r],r]r]'j2r] 


2++ 


4* 


1934 ± 20 


271 lb 25 


0.0 zb 0.08 


1462 


7r7r,r7r?,/2??,a27r,GAMS,VES 


2++ 


3* 


2010 ± 25 


495 ± 35 


1.51 ± 0.09 


694 


rjr] X 2,WA102 


2-+ 


4* 


2267 ± 14 


290 ± 50 




1349 


f2V',[f2V]L=2,aoT^, /o(1500)r7 


2-+ 


3* 


2030 ± 16 


205 lb 18 






[a27r]L=o,2,ao7i" 


2-+ 


4* 


1870 ± 16 


250 ± 35 






/277,a27r,ao7r,WA102 


1++ 


2* 


2310 ± 60 


255 ± 70 




882 


f2V, f2V', crri 


1++ 


3* 


1971 ± 15 


240 lb 45 




1451 


/2?7, a27r aovr, {fo{980)r],ar],ar]') 


0++ 


3* 


2337 ± 14 


217 lb 33 






vr 71,7777 


0++ 


4* 


2102 lb 13 


211 lb 29 






777? X 2,777?',Mark III,BES I 


0++ 


3* 


2020 ± 38 


405 lb 40 






TTTT, WA102(7r7r,47r) 


0-+ 


3* 


2285 ± 20 


325 ± 30 




1342 


/o(1500)77 X 2,{av),iav') 


0-+ 


2* 


2010li 


270 ± 60 




1189 


a27r,/2??,o'77, BES II 


Table 2 



/ = 0, C = +1 resonances . Errors cover systematic variations in a variety of fits 
with different ingredients. The sixth column shows changes in S" = log likelihood 
when each component is omitted in turn from the fit to rjir^ir^ data and remaining 
components are re-optimised. In the last column, channels in parentheses are weak 
compared with the others, or less reliable. 



Note the difference in sign in the last two of these expressions for T-matrix 
elements Tj+i^j. Differential cross sections and polarisations are given by 

AoNda/dn = Tr{F*(7YF) = 27m(F;+F+_) (9) 
(i(7/dQ=|F++|2 + (10) 
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Fig. 11. Differential cross sections for pp from 360 to 1300 MeV/c, com- 

pared with the fit; panels are labelled H for data of Hasan et al. and E for those of 
Eisenhandler et al. Each distribution is normalised so as to integrate to 2it. 

The polarisation is very sensitive to the difference between amplitudes for 
L = J ± 1, for example and amplitudes. It is also phase sensitive. 

The analysis of 2-body data alone [13] leads directly to a unique and accurate 
solution for all 0+, 2+ and 4+ states. Intensities of fitted partial waves are 
shown in Fig. 14. A small detail is that the low-mass tail of a higher 6+ 
resonance is included, though the mass of that state is above the available 
mass range and is therefore not well determined. The /4(2050) is observed 
strongly in Fig. 14 in all of tt'^tt^, rjrj and 77?]', as well as in ryvr^TT^ data discussed 
below. The 2-body data required the presence of two states, two and 
two 0+ states. The strong peaking of the ^P2 amplitude at the lowest mass is 
associated with /2(1934), or /2(1910) of the PDG. Fig. 9. of Ref. [13] shows the 
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Fig. 12. As Fig. 11, 1351 to 2230 MeV/c. 

effect of dropping each of tliese states one by one from tlie analysis. Removing 
any of tliem causes clear visual disagreements with the data as well as large 
changes in log hkelihood. 

A cross-check on the analysis is possible between data and vr vr"'". States 
with / = 1, C = — 1 contribute to vr"?!"*" but are determined separately from 
LJTi data discussed below. If one allows for their small intensity in n'n^ and 
folds the angular distribution forward-backwards about 90°, in order to remove 
interferences between / = and 1 = 1, the resulting angular distribution 
agrees accurately in shape with and their absolute magnitudes agree 

within 1.1±2.3%. This is a valuable check on the entire process of selecting and 
reconstructing Crystal Barrel data, and the correction for rate dependence. 
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Fig. 13. Polarisation data for pp ^ tt 7r+, compared with the fit. 
5.2 3-body data 

Resonances with = 0~, 1+, 2^, 3+ and 4^ appear only in 3-body channels. 
However, statistics for 777r°7r° are very high, typically > 70 K events per mo- 
mentum. There are two dominant channels, f2{1270)r) and a2(1320)7r which 
are both well determined. Mass projections at some momenta are shown in 
Fig. 15. Small /o(980), /o(1500) and ao(980) signals are also visible. 

The /4(2050) acts as an interferometer at the lower momenta, and hkewise 
/4(2300) at higher momenta. Phases of many other partial waves show rather 
small mass dependence with respect to these f^. This is illustrated in Fig. 16; 
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Fig. 14. Intensities for partial waves for tt^tt^, r)r) and rjrj'. Dashed curves show 
contributions from L = J — 1 and dotted curves from L = J + 1. 

there, M refer to the heUcity in the initial pp state. It is unavoidable that reso- 
nances need to be introduced into S"*", l"*" and 0~ partial waves. The behaviour 
of 2+ partial waves is somewhat more complex. The 2-body data define very 
precisely the existence of states which are dominantly at 1934 and 2240 
MeV and at 2001 and 2293 MeV. The first of these may be identified with 
the /2(1910) observed by GAMS [21] and VES [22] collaborations. The small 
spacings of 2^ states leads to rather rapid variations of phases with mass, but 
the 3-body data are entirely consistent with 2-body results. 

There are then two lucky breaks. The first is that data for the final state 
ryV^TT*' are dominated by a single resonance in /2(1270)?7' [23]. Dalitz plots are 
shown in Fig. 17. There is an obvious band due to /2(1270) along the lower left- 
hand edge of these Dahtz plots. A very clear peak is observed in the integrated 
cross section for f2r}' at 2267 MeV, see Fig. 18. It is consistent with a resonance 
required near this mass in rjir^n^, but gives a better definition of resonance 
parameters. In turn, this 2~ state then acts as a powerful interferometer to 
determine other partial waves at high masses. 
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Fig. 15. Mass projections on to M(7r7r) and M{'Krj) for rj-K^ir^ data at (a) and (b) 
900 MeV/c, (c) and (d) 1200 MeV/c, (e) and (f) 1525, (g) and (h) 1800 MeV/c; 
histograms show the result of the fit. Data are uncorrected for (small) variations of 
acceptance; these are included in the maximum likelihood fit. 

The second lucky break comes from 777777 data [24] . Although statistics are low, 
there is a strong peak at high mass due to the well identified /o(1500)77 S-wave 
final state with quantum numbers 0~. This is shown by the curve labelled 0~ 
in Fig. 19; in plotting this curve, interference between the three /o(1500) bands 
has been removed. There is a small but distinctive contribution from the 2~ 
state at 2267 MeV in the /o(1500)77 D-wave. The analysis is very simple and 
leads to a cleanly identified 0~ resonance at 2285 MeV. 

Returning to rjmr data, there remain a 4~ state, two 1^ states and a lower 0~ 
state which are needed in the full combined fit. 
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2* Uri with L=1, M=0 2* o^n with L=1, M = 4* ^ f^Tj with L=3, M=0 
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Fig. 16. Relative phases of partial waves for -qir^iT^. The phases for 2+ and 4+ with 

M=0 are relative to 4^ a27r with L=3, M=0; the phases for l"*", S"*" and 4"*" with 
M=l are relative to 4+ — > /2?7 with L=3, M=l; the phases for 0~ and 2~ are relative 
to 2~ /2?? with L=0. 



5.5 2" spates 



It is necessary to discuss separately 2~ states at 1645, 1860 and 2030 MeV. 
They appear in t^tt^tt^tt^ data in the production reaction pp TT^{ri7r7r). This 
was first studied [25] in data at just two momenta: 1200 and 1940 MeV/c, prior 
to the running of the other momenta . Two 2~ / = states were observed. The 
lower one was initially observed decaying purely to a2(1320)7r at 1645 MeV. 
This fits nicely as the expected partner of the well known 1 = 1 7r2(1670). 
However, the same data contained a second strong signal in the /2(1270)?7 
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Fig. 17. Dalitz plots at all beam momenta for ry'-zr^Tr^ data. Numerical values indicate 
beam momenta in MeV/c. 




Q ^ 

1.9 2 2.1 2.2 2.3 2.4 

y (GeV) 



Fig. 18. The integrated cross sections for the final state rj'TT^ir^, after correction for 
backgrounds and for all decay modes of ry', rj and ir^; the full curve shows the fit 
from the amplitude analysis; the dotted curve shows the arj' contribution. 

final state at 1870 MeV. The second signal could not be explained as the high 
mass tail of 772(1645) decaying to f2r]. Furthermore, the relative strengths of 
the two signals were quite different in data at 1200 and 1940 MeV/c. Two 
separate states were required. 

Both resonances were subsequently confirmed by the WA102 collaboration. 
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M(VVV) (GeV) 

Fig. 19. The fit to the mean integrated cross section for pp 3r] from 67 and 
IO7 data, (b) r/(2285) ^ [/o(1500)r/]L=o (dashed), and ?72(2267) ^ [/o(1500)7?]l=2 
(chain curve). 

They observed the 772 (1645) in decays to rjTnr and 772 (1870) in decays to both 
rjUTT and KKn [26]. However, the KKn decay mode is not strong enough 
to make an interpretation hkely as the ss partner of 772 (1645). A feature of 
pp annihilation is that well known ss states such as /2(1525) are produced 
very weakly, if at all; the 7^2(1870) is produced rather strongly. A possible 
explanation is that it is the radial excitation of 7^2(1645), but that would be 
expected to he around 2000-2050 MeV, near /4(2050). 

In later data, a much more extensive study was made of 777r°7r°7r° data with 
statistics a factor 7 larger [27]. A further 2~ state appears as expected at 2030 
MeV, decaying dominantly to the a2(1320)7r D-wave. 

The main features of the data are shown in Fig. 20 as scatter plots of M(7r7r) v. 
M^nr]) for three rjTnr mass ranges centred on 7^2(1645), 7^2(1870) and 7/2(2030). 
In (c), there is a vertical band due to 02(1320) from the decay of 7^2(1645); 
the horizontal band at 1000 MeV is due to the low mass tail of the stronger 
772(1870) — > f2r]. In (b), the horizontal band from 7/2(1870) is dominant and the 




Fig. 20. Scatter plots of M{Tnr) v. M{r]'iT) for three ranges of M{r]inT): (a) 1945-2115 
MeV, centred on r/2(2030), (b) 1775-1945 MeV over r?2(1870) and (c) 1560-1750 
MeV, centred on 7/2(1645). 
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Fig. 21. S = -log likelihood v. masses and widths of (a) and (b) 7/2(1645), (c) and 
(d) r?2(1875), (e) and (f) 772(2030). 

vertical a2(1320) band becomes weaker. In (a), the /2(1270) becomes weaker 
and there is a strong peak at the intersection of /2(1270) and a2(1320); the am- 
plitude analysis shows that the 7^2(2030) decays dominantly to the a2(1320)7r 
D-wave. The essential point of the analysis is shown by the branching ratios 
in Table 3. [These refer directly to data and are not corrected for other charge 
states or for other decay modes of the 02 and 77.] Branching ratios change 
dramatically with mass, requiring three separate states. The /2(1870) decays 
dominantly to f2i] and aQir; the 7/2(2030) has dominant 02(1320)7/ decays with 
L = and 2, and the f2r] decay is weak. 

Fig. 21 shows log hkelihood for the three 2~ states for mass and width. There 
is some evidence for decay to rja, though that is a difficult final state to identify 
because the a is broad. 

The rj2 states couple to pp ^D2. Prom the straight-hne trajectories of Fig. 1, 
it appears that 7/2(1870) is an 'extra' state. Isgur, Kokoski and Paton [28] 
predicted the existence of hybrids around 1800-1900 MeV. 

The hybrid consists of q and q connected by a gluon having one unit of orbital 



Resonance 


7/2(1645) 


%(1870) 


7/2(2030) 


BR[ao'!r]/BR[a2Tr] 


0.21 ±0.13 


3.9 ±0.6 


0.21 ±0.02 


BR[f2r]]/BR[a2Tr] 




4.5 ±0.6 


0.25 ± 0.09 



Table 3 

Decay branching ratios for 7/2(1645), 7/2(1870) and 7/2(2030). 
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angular momentum about the axis joining the quarks. It rotates hke a skipping 
rope between q and q. It is predicted that this orbital angular momentum will 
be transferred to one pair of quarks in the final state [29,30]. The expected 
decay modes are therefore to one combination of quarks and one S-state 
combination. Favoured decay modes are /2(1270)?7 and a2(1320)7r, in agree- 
ment with 7^2(1870) decays. However, the f2r] decay mode is strikingly strong. 
One can use SU(3) to calculate the expected ratio BR[a2Ti]L=o/ BR[f2ri\, in- 
tegrating their available phase space over 02 and /2 line-shapes. If it is a qq 
state, the predicted ratio for ?72(2030) is 9.7 after correcting for all charged 
states and decay modes of 02 and 77; this compares well with the observed 
10.0 ± 1.8. However, for 7^2(1870), the prediction is 8.5 compared with the ob- 
served 1.27 ±0.7. The conclusion is that the decay to /2?7 is distinctive and far 
above the value expected for a qq state. It will be shown below that there is 
an / = 1 partner at 1880 MeV decaying to a2(1320)?7. If these make a hybrid 
pair, it is the first pair to be identified. Further hybrids containing hidden 
strangeness are then to be expected at ~ 2000 MeV decaying to ^2(1430)77 
and at ~ 2120 MeV decaying to /^(1525)77 and/or X2(1430)X. 

5.4 Star ratings of states 

All states of Table 2 have a statistical significance > 25 standard deviations. 
The even spin states are even more significant in 2-body decays. Argand dia- 
grams for 777777 are shown in Fig. 22 and for 7r7r decays in Fig. 23. 

The data require a larger and more linear phase variation than can be produced 
by a single resonance. In electronics, it is well known that a linear phase 
variation with frequency may be obtained with a Bessel filter [31], which has 
a sequence of poles almost linearly spaced in frequency. Physically, such a 
system generates a time delay which is independent of frequency, since group 
velocity is proportional to the gradient of phase v. frequency. It appears that 
a sequence of meson resonances likewise conspires to produce approximately 
a linear phase variation with mass squared s. 

Intensities of the main individual partial waves in t^tt^tt^ are displayed in Fig. 
24. In each panel, titles indicate in descending order full, dashed, chain and 
dotted curves. 

The /4(2050) is seen as a clear peak in tttt, 7777, and a2(1320)7r, Fig. 24(g) and 
is clearly 4*. It peaks at 2070 MeV because of the effect of the centrifugal 
barrier for production from pp, but the mass of the Breit-Wigner resonance is 
substantially lower, at 2018 MeV. The state at 2283 MeV is the /4(2300) of 
the PDG. It is not conspicous as a peak in 2-body data, but has been found by 
several earlier partial wave analysis of 7r~7r"'" data by Hasan [1], A.D. Martin 
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Fig. 22. Argand diagrams for prominent partial waves in -qn tt data; crosses mark 
beam momenta; the amplitudes move anti-clockwise with increasing beam momen- 
tum; in (a) and (b) dashed curves show the effect of removing the lower 0~ resonance. 

and Pennington, [32], B.R. Martin and Morgan [33] and Carter [34]. It need 
not appear as a peak in pp — > tttt because of interference with /4(2050). It does 
however appears clearly as a peak in f2i] in Fig. 24(g) at 2320±30 MeV. Recent 
VES data also display conspicuous 4"*" peaks in lulu at 2325 ±15 MeV with 
r = 235 ± 40 MeV [35] and in rjirn data with M = 2330 ± 10{stat) ± 20{syst) 
MeV with T = 235 ±20 ±40 MeV [36]. Peak positions agree well with Crystal 
Barrel. This is also a 4* resonance. 

The 4~ state {pp ^^4) appears only in rjim data. Amphtudes are small, but 
rather stable in all fits. The dominant decay is to ao(980)7r. Fig. 24(m). This 
resonance is classified as 2*. Log likelihood changes by 558 when it is omitted. 
General experience is that any contribution improving log likelihood by more 
than 500 has well determined parameters; any affecting log likelihood by > 300 
is definitely needed. 
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Fig. 23. Argand diagrams for 7r"'"7r amplitudes; crosses mark masses at 100 MeV 
intervals of mass beginning at 1900 MeV. 

The upper 3+ state decays to /27r with both L = 1 and 3; it produces the 
strongest intensity at high mass in rjinr data, Fig. 24(d). It also appears as a 
clear peak in rj'n^n^ data. The Argand diagram of Fig. 22(h) shows two loops 
requiring two 3^ resonances. The lower one appears clearly as peaks in Fig. 
24(d) and (1) in /2(1270)r7 and ao(980)7r and more weakly in a2{1320)n . Phases 
are well determined by interference with /4(2050), so it really should be given 
4* status. Its mass is well determined with an error of only ±8 MeV. 

The upper 2^ state at 2267 McV appears as a clear peak in the raw rj'Tnr 
data of Fig. 18. It is also observed clearly in rjirn in decays to [f2i]]L=2 in Fig. 
24(c)(chain curve) and to both ao(980)7r and /o(1500)?7, Fig. 24(k). It deserves 
to be classified as 4*. The Argand diagram of Fig. 22(e) requires a lower 2~ 
state also. It appears as a peak in Fig. 24(c) in both [f2v]L=o and a27r; it is 
also seen in several decay modes in ?77r°7r°7r° data. So its existence is not in 
doubt. Because of its overlap with 7^2(1870), its mass and width may have 
some systematic uncertainty, so it falls into the 3* category. 

The first four 2+ states of Table 2 are all extremely well determined in 2-body 



33 



20 
15 
10 
5 



0^ 
1.9 



^ 

'So-f,(1 275)7) L = 2 
'S[,-Oj(1320)tiL = 2 



_-fj(l 275)7] L=0 
l2-aj(1320)nL=0 40 
l2-f2{1275)7| L=2 

1320;nL=2 30 





-fj(1275)77 L=1 




-a,(1320)7iL=1 




-■'2(1275)7; L=3 


- T., 


-ajC32CkiL=3 




x2 / \^ 




/ Y'- . X • 




/ X ■ x2 


: / 

- 





2.17 2.44 1.9 



60 
40 
20 



Pz-fi(1275)7i L = 
^P2-a2(1 320)7iL = 
^P2-f2(1275)7i L = 3 
-\P2-<i!(1320)7iL=3 




S„-=(1G00)t, l = q 
' S„- a„f98G)TT L = Q 
'S„- ft;975)-, L = Q 
's„-f.{;5C0)?i L = u 





2.44 1.9 



30- 



- 'p, 


-tj(IOOO)i) L=1 


- 'p, 


- Oo{980)7T L=l 


- 'p, 


- f,(975)7) L=l 


- 'p 


-o,(1450)tiL=1 







2.44 1.9 



2.17 



: 'p2 




1 640' 






-□2 


1 640) 






-02 


1 O-IO^ 


77 1=3 


; C2 


-□2 


1 640 


77L=G 








, ^2 






' 










v2 






)\/Y 





2.44 1.9 



2.17 2.44 1.9 

M(pp) (GeV) 



2.17 



2.44 1.9 

1.5 
1 

0.5 


2.44 1.9 



2.17 2.44 

(p) 



p, -02(1640)711= 

F3-a2(' 640)77 1= 
- 'H4-02(' 610)771= 



2.17 



2.44 



Fig. 24. Intensities fitted to individual partial waves, and corrected for all r] and 
7r° decays. The titles in each panel show in descending order full, dashed, chain 
and dotted curves; (7(1000) denotes the broad /o(400 — 1200) in the tttt amplitude. 
Numerical values indicate factors by which small intensities have been scaled so as 
to make them visible. 

data and r/)/r7r and are all 4*. The /2(1980) will be discussed later as /2(1950) 
of the PDG. 

Coming to low spins, states are less well defined. The lower 1"*" state i^Pi) 
is large in several channels in Figs. 24(b) and (j), dominantly /2?7, a2Ti and 
ao(980). Despite the fact that it lies at the bottom of the mass range, parame- 
ters are quite stable: M = 1971 ± 15 MeV, F = 241 ±45 MeV. It just deserves 
3* status. However, the upper 1+ state is the least well determined state in rjnn 
data. The reason is that it makes only small contributions to the distinctive 
channels 02(1320):^- and /2(1270)?7, and there arc large interferences with other 
amplitudes. As the resonance mass is changed, these interferences are able to 
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absorb the variation with small changes in log hkelihood. Nonetheless it does 
improve log likelihood by 882, a very large amount. Since it also appears in 
Tj'TTTT, it qualifies as 2*. 

The upper 0~ state {pp ^Sq) at 2285 MeV appears as a clear peak in the raw 
3r) data of Fig. 19. It appears again in the /o(1500)?7 channel in rjinr with a 
magnitude consistent with the known branching ratio of /o(1500) between tttt 
and rjT]. It also contributes to arj. It is classified as 3*. The lower 0~ state 
is close to the pp threshold and its amphtude goes to a scattering length 
at threshold. This makes it difficult to separate from a threshold effect. On 
the basis of Crystal Barrel data it is only 1*. However, new BES II data on 
J/^ — > 747r identify this state clearly at 2040 MeV in An [37], raising its status 
to 2*. 

Finally we come to O"*" states. The upper 0^ state at 2337 MeV is remarkably 
well defined in all of tt~71~^, tt^tt^ and rir]. This qualifies it as 3*. The /o(2100) 
appears as an rji] peak in t^t^tt data at 1940 McV/c [38], see Fig. 14 of Section 
5.1. It is also conspicuous in rjrj as shown in Fig. 14 and contributes to rjr]'. 
If it is dropped from the 2-body analysis, increases by 4030, an enormous 
amount. Allowing for its low multiplicity (2J+1), it is one of the most strongly 
produced states in pp annihilation. It was discovered in J/\E' — s> '-fAir [39] and 
also appears as a very clear peak in E760 data on pp 77777?" [40], but has no 
spin determination there; these data are illustrated below in Fig. 25. Its mass 
and width are very well determined by the E760 data. Crystal Barrel data 
agree within a few MeV with that determination. It is classified as 4*. 

An /o(2020) with a width of 400 MeV has been reported by the WA102 col- 
laboration [41] in decays to both 27r and An. Adding this to the analysis of 
Crystal Barrel 7r°7r° data gives a significant improvement in of 370 and a 
visible improvement in fits to some da/dQ data. However, overlap and inter- 
ference with the strong /o(2100) make a determination of its mass and width 
difficult. It also obscures /o(2100) in nn in Fig. 14. From Crystal Barrel data 
alone it is a 1* resonance, but its observation in two channels of WA102 data 
raises its status to 3*. 

In summary, all states of Table 2 are securely determined in mass except for 
the upper 1+ state at 2310 MeV and the lower 0" state at 2010 MeV; both 
are certainly required, but have rather large uncertainties in mass, ±60 MeV. 

5.5 /o(1710) and /o(1770) 

Data from MARK HI [42] and DM2 [43] on J/^ -fKK establish the ex- 
istence of /o(1710) decaying to KK It is also observed clearly in data on 
J/^ — > ujKK [44-46] . From these and the absence of evidence for decays to nn, 
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one can deduce a branching ratio [/o(1710) — > 7r7r]/[/o(1710) KK] < 0.11 
with 95% confidence. The /o(1710) is an obvious candidate for an ss state 
completing a nonet with /o(1370), ao(1450) and i^ro(1430). 




(tt'tt-ttV) mass (GeV/c') MASS (?] 77) GeV/c' 

Fig. 25. (a) Mark III data for J/^ — > jAtt, (b) the rjr] mass spectum from E760 data 
on pp rjijiT^ at 3.0 GeV. 

A radial excitation of /o(1370) is expected in the same general mass range. 
There is ample evidence for this state, and the PDG used to call it X(1740). 
Today, all entries are lumped under /o(1710). However, these need to be ex- 
amined critically. In an analysis of Mark III data for J/^^ — > ■jTr~^Tr~TT~^7i~ , 
Bugg et al. identified a clear peak in An at 1750 MeV as a 0"*" resonance [39]. 
The Mark III data are shown in Fig. 25. If it is identified as /o(1710) it would 
require a branching ratio for J/\E' 'yAn of (10.9 ± 1.3) x 10~^, larger than 
that to 'jKK, (8.5tJ;9) x 10"''. That does not fit nicely with an ss state, and 
suggests the presence of a second state. 

In J/^f — > ^TTTT there is a peak at 1765 MeV but no corresponding signal in 

J/^ (j)KK [46]. New data from BES confirm and strengthen this observa- 
tion [47]. These data arc shown in Figs. 26 and 27. There is a clear peak in 
Fig. 26(c) and (e) at 1780 MeV; = 0"*" is strongly favoured over 2+. In Fig. 
27 there is no corresponding peak in (f)KK, but a small shoulder in Fig. 27(c) 
at 1650 MeV; Mark III results are similar. This small shoulder is fitted by 
/o(1710) interfering with /o(1500) and /2(1525). The branching ratio nir/KK 
for the 1765 MeV peak in ^tttt is at least a factor 25 larger than for /o(1710), 
so two separate states are definitely required. The natural interpretation of 
the peak in ^tttt is in terms of a second /o(1770) decaying dominantly into tttt 
and An but not significantly into KK. This fits in with the second 0+ state 
expected in this mass range. There is however a long-standing puzzle why the 
/o(1710), which seems to be dominantly ss, should be produced in uKK while 
the /o(1770), apparently mostly non-strange, should be produced in the final 
state (j)KK. 

In Crystal Barrel data on pp rjrjn^, there is evidence for the same /o(1770) 
decaying to rji] [48]. Fig. 28 shows the rjr] mass projection for data at 900, 
1050 and 1200 MeV/c. The histogram shows the best fit without inclusion 
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Fig. 26. BES II data from J/'^ — > (pTT^7r~ . (a) and (b) measured and fitted Dalitz 
plots; (c) and (d) tttt and 4>tt mass projections; the upper histogram shows the fit 
to data and the lower shaded histogram show background derived from a sidebin to 
the (p; in (d) the dip at 1.5 GeV arises from a cut to remove K*K*; (e) and (f) fitted 
O"*" and 2"*" contributions; the dashed histogram in (e) shows the a contribution. 

of /o(1770). At all three momenta there is a small surplus of events, which 
may be fitted well by including a further contribution from /o(1770) with 
M = 1770 ± 12 MeV, T = 220 ± 40 MeV; these values are in good agreement 
with data discussed above. The first three panels of Fig. 29 show log likelihood 
V mass at the three momenta. Full curves are for the spin assignment and 
dashed curves for spin 2; spin is obviously preferred. The last panel shows 
log likelihood v. width. 

The /o(1710) could decay to rjr] via their ss components. However, ss final 
states are conspicuous by their absence in pp annihilation. In the same rirjir^ 
data, a scan of the 2"*" amplitude v. mass shows no peaking at /2(1525), so 
/o(1710) appears an unlikely source of the rjr] signal. Furthermore, mass and 
width fitted to /o(1770) appears significantly different to those of /o(1710). 

A somewhat different fit to these data has been reported by Amsler et al. 
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Fig. 27. As Fig. 28 for J/^- (pK+K'. 




Fig. 28. Data for pp — > rjrjir^. Histograms show the fit without /o(1770) at (a) 900 
MeV/c, (b) 1050 MeV/c and (c) 1200 MeV/c. 

[49]. This fit omits /o(1770) and generates the 1770 MeV shoulder as the 
tail of an interference between /o(1370) and /o(1500); since they lie close to 
one another, they behave like a dipole with a long tail at high masses. The 
difficulty with this fit is that the /o(1370) branching ratio to 7777 is known 
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to be small compared to that for nn; the large /o(1370) signal required is 
considerably larger than is fitted to /o(1370) in Stt data at all momenta, see 
Section 6 below. 




Fig. 29. Log likelihood v. mass of /o(1770) full lines J = 0, dashed lines J = 2, (a) 
at 900 MeV/c, (b) 1050 McV/c and (c) 1200 McV/c; (d) log likelihhod v. width at 
1050 MeV/c (full curve) and (d) 1200 McV/c (dashed). 

The first evidence for a signal in rir] at 1744 ±15 MeV was reported by the 
GAMS collaboration [50,51]. It was observed to have an isotropic decay angu- 
lar distribution, consistent with = 0~^. 



5.6 SU(3) relations for 2-body channels 



Amplitudes for pp — > n^n^, rjrj and rjrj' throw valuable light on the quark 
content of resonances. It is well known that rj and rj' may be expressed as 
linear combinations of \uu + dd > /\/2 and ss. Denoting the former by nn, 
their quark content may be written in terms of the pseudoscalar mixing angle 
© as 

|nn > = cos 0|?7 > + sin ©I?]' > (11) 
|ss > = — sin ©|?7 > + cos ©I?]' >, (12) 

where cos© ~ 0.8 and sin© ~ 0.6. Amplitudes for decay to tt^tt^, rjr) and 7777' 
may then be written compactly as [13] 

Ivr'^Tr'^ > cos^ © 

f{nn) — p 1 ^^[7777 >+ cos ©sin ©I 7777' >, (13) 

V 2 V 2 
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/(ss) = V^[sm^© I r^r^ > +cos^ Q\r}r] > -v^cos © sin ©jr^r^' >]. (14) 

The second of these equations includes a factor A. This factor has been fitted 
empirically by Anisovich to a large variety of data on strange and non-strange 
final states at high energies [52] . Its value is y/X — 0.8-0.9 and the central value 
0.85 is adopted in the Crystal Barrel analysis. 

Resonances R observed in data may be linear combinations of nn and ss: 

R — cos ^\qq >+ sin ^\ss > . (15) 
Amphtudes for the decay of R to these three channels are given by 

/(A°) = ^^ (16) 

COS ^ 

/(7;7;) = ^_(cos^0 + V2Asin2 0tan$) (17) 
v2 

/(r;r7') = ^^^cos©sin©(l- V2Atan$). (18) 
v2 



Most flavour mixing angles $ optimise close to zero, i.e. close to pure nn 
states. Some mixing is essential, since r/r/ differential cross sections are visibly 
different to tt'^tt'^. The small mixing angles allow modest changes in branching 
ratios between 7r°7r°, rjr) and rjr)'. However, there is a striking exception for 
/o(2100). This requires $ = (65 ± 6)°. The reason for this result is that the 
ratio of observed amplitudes of /o(2100) is 

ttV : rjT] : rjrj' = 0.71 ± 0.17 : 1 : -0.85 ± 0.45. (19) 

For an unmixed nn state the branching ratio BR[7r^7r^]/ BR[r]r]] would be 
0.8~^ = 2.44. The /o(2100) is one of the strongest 0+ resonances. It makes up 
(38 ± 5)% of the rjr) signal. Intensities of rjr) and rjr)' cross sections v. mass are 
shown in Fig. 10(c) and (d). The /o(2100) has the anomalous property that it 
is produced very strongly in pp yet it decays dominantly to ss. Below it will 
be considered as a candidate for a 0"*" glueball predicted close to this mass. A 
glueball would have a mixing angle of +35.6°. The departure of the fiavour 
mixing angle $ from this value can be accomodated by some mixing with a 
nearby ss state. An overall picture of the fit to 7r°7r°, rjr) and rjr]' data in shown 
in Fig. 30. 



6 / = 1,C = +1 resonances 



The channels available are tt^tj, tt^tj', 37r° and rjrjn^. Statistics are very high 
for 37r°, typically > 250 K events per momentum; backgrounds are very low 
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Fig. 30. Differential cross sections for pp — 7r°7r° (first two rows), r/r/ (rows 3 and 
4) and 7777' (bottom 2 rows), compared with the combined fit (full curves). For 7777', 
the average of 47 and 87 data is used. 

[53]. Since u and d quarks are nearly massless, one expects the spectrum of 
/ = 1, C = +1 mesons to be close to that for / = 0, C = +1. Although results 
are indeed similar, they are very much less accurate. Some resonances are well 
defined, but several are not. The well defined resonances agree reasonably well 
with / = 0, C = +1. The remainder fall into line with 7 = if one chooses 
one of two alternative solutions to the data. 

The fundamental difficulty is that there are two discrete solutions for n^rj and 
7r°?7' [12]. To get from one solution to the other, two 2+ amplitudes must be 
flipped in sign; a substantial change also appears in ^7^4. The reason for the 
ambiguity is the lack of polarisation data, which are needed to separate 
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and cleanly; the phase sensitivity of the polarisation data is also missing. 

The 37r° channel is dominated by /2(1270)7r. Two 4+ resonances appear reason- 
ably clearly as does a particularly strong 3+ resonance at 2031 MeV. Further 
analysis requires resonances in all partial waves up to 4~, in close analogy to 
I — 0, C — +1. However, determinations of resonance parameters are rather 
fuzzy. An initial analysis was made of the /2(1270)7r channel alone [53]. 

This analysis was made at single energies, and resonance parameters were 
deduced by fitting curves to magnitudes and phases. This technique is haz- 
ardous, since magnitudes and phases contain non-analytic noise at individual 
momenta. Where relative phases between partial waves arc close to or 180°, 
phases can depart seriously from physical values and can distort conclusions. 
Removing this noise is important. The later analysis of Ref. [4] was made to 
all momenta simultaneously, fitting resonances directly to data. There were 
substantial changes to the solution, particularly at 600 and 900 MeV/c. 

There are several reasons for this instability: 

- (i) As well as /27r in Btt" data, there is an almost flat, featureless physics 
background. It may be fitted with the a pole at 550 MeV and a broad, slowly 
varying tttt amphtude over the mass range 1-2 GeV, well known in tttt elas- 
tic scattering [16,54]. The entire tttt S-wave amphtude will be denoted by a 
as a shorthand. The cttt amplitudes appear mostly in pp singlet states. They 
interfere with /27r singlet amplitudes and the interference leads to a poor 
definition of resonance parameters; substantial shifts of masses and widths 
may be accomodated by changes to phase angles. 

- (ii) Lack of polarisation data leads to uncertainty in r2 and parameters, 
i.e. separation between and and between ^F^ and ^H^. 

- (iii) There is a rapid change in angular distributions between the two lowest 
beam momenta: 600 MeV/c (M = 1962 MeV) and 900 MeV/c (M = 2049 
MeV) ; this is shown on Fig. 2 of Ref. [4] . The experiment was supposed to 
run also at 450 and 750 McV/c, but these momenta were lost because of the 
closure of LEAR. For / = 0, the data from PS172 extend the range down 
to 360 MeV/c (M = 1912 MeV) in steps of 100 MeV/c. 

The ultimate conclusion is that 37r^ data do not define the partial waves 
cleanly, despite the high statistics. 

Amongst the 3-body data are typically 5000-9000 rjrjTi^ events per momentum. 
It turns out that these data arc actually more valuable than 37r*^. In rirjiT^, there 
are strong contributions from spinless channels ao(980)?7, /o(1500)7r, /o(1770)7r 
and /o(2100)7r. For these channels, amplitudes depend only on Legendre poly- 
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nomials for the production angle. The resulting analysis is very simple. Singlet 

states dominate. As a result, 0~, 1+, 2~, S"*" and 4~ resonances emerge well 
defined. There is also a further distinctive state at 2255 MeV. This infor- 
mation feeds back to stabilise the fit to Stt^. 

The final element in the analysis is the observation of a 7r2(1880) state, — 
2~, in the production process pp 7r2(1880)7r, tt2 tji^tt in rjrjTi^Ti^ data. 
This state overlaps partially with a further 2^ state at 2005 MeV. They can 
actually be separated by the fact that 7r2(2005) decays to ao(980)r; whereas 
the 7r2(1880) does not. However, for other decay channels there is potential 
ambiguity between these two close resonances. 
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Fig. 31. Angular distributions for w r] corrected for acceptance. Black circles show 
results from 47 and open triangles results from 87. Full and dashed curves show fits 
from solutions 1 and 2 respectively. Dotted curves illustrate the acceptance. Open 
squares show earlier results of Dulude et al at neighbouring momenta, after renor- 
malisation to obtain the optimum agreement with Crystal Barrel normalisation. 

6.1 2-body data 



Angular distributions for Trrj are shown in Fig. 31 [12]. There is again excellent 
agreement between 47 data (black circles) and 87 (open triangles) in both 
shapes and absolute normalisations. Open squares compare with earlier data 
of Dulude et al. [55] after renormalising them to agree with the normalisation of 
Crystal Barrel data. Angular distributions for 71^77' are shown in Fig. 32. There 
are typically 1000 events per momentum. Here one sees differences between two 
alternative solutions, shown by full and dashed curves. Intensities of partial 
waves from these two solutions are shown in Fig. 33. The main difference lies 
in 2+. 

Both solutions require the resonances at 2255 and 2030 MeV and some 



44 




0.5 10 0.5 10 0.5 1 




0.5 10 0.5 1 0.5 

CosO CosO CosO 



Fig. 32. Angular distributions for ir'^r]' corrected for acceptance. Full and dashed 
curves show solutions 1 and 2 respectively. Dotted curves illustrate the acceptance. 

contribution. In the earliest analysis of Ref. [12] there was no input from 
37r° data and the upper resonance at 2175 MeV was not included. Even 
in the final combined analysis, this is the least significant contribution. 

Both solutions require a 0+ state near 2025 MeV, though its mass and width 
have sizable errors. A lucky break is that independent evidence for this state 
appears in rjir^n'^Tr'^ data. A small but definite ry(1440) signal is observed there 
in rjTrn [56]. It is illustrated in Fig. 34. This channel may be fitted to production 
amplitudes using simple Legendre polynomial for L = 0, 2 and 4. The result 
is complete dominance by L = 0, with no significant evidence for any L — 2 
or 4 amplitudes. This requires production from an initial I — 1, J^^ — 
initial state. The cross section for production shows a distinct peak at 2025 
MeV, Fig. 35. It may be fitted with a resonance at 2025 MeV, F = 330 MeV 
after inclusion of the pp L = 1 centrifugal barrier for production. 

However, even after this lucky break, the 2-body data fail to locate the further 
state expected in the mass range 2300-2400 MeV as partner to /o(2337). Note 
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Fig. 33. Intensities of 0+, 2"*", 4+ and G"*" partial waves for pp ir^r] for the better 
solution (upper row) and poorer solution (lower row). Dashed curves show intensities 
for and ^^4, dotted curves for and ^i?4, full curves their sum. 

that the oq expected between 1450 and 2025 MeV is also missing. As a result, 
the 0"*" sector for 7 = 1 is poorly known. 
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Fig. 34. Mass projections for all r]-K^TT^ combinations in rjSir^ data after selecting 
events where M(r/7r) lies within ±50 MeV of ao(980) or alternatively M^(7r7r) > 0.6 
GeV^. Histograms show fits from the amplitude analysis. 
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Fig. 35. Cross sections for pp ry(1440)7r with L = 0. The curve shows the cross 
section for a resonance with M = 2025 MeV, F = 330 MeV, including the pp i = 1 
centrifugal barrier at threshold. 
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Fig. 36. Dalitz plots for pp ryryvr; numbers indicate beam momenta in MeV/c. 
6.2 pp ^ TjfjT:'^ 



These data are discussed in detail in Ref. [5]. Dalitz plots are shown in Fig. 36 
at all momenta. There are obvious vertical and horizontal bands at M'^{ri7i) = 
1 and 1.7 GeV^ due to ao(980) and a2(1320). The is also a diagonal band at 
the highest four momenta due to /o(1500). 

The TTTT mass projection is shown in Fig. 37. The strong /o(1500) is clearly 
visible. There is also a shoulder due to /o(1770) at beam momenta 1200-1642 
MeV/c. There is one other important feature of Fig. 36. At the intersection 
of vertical and horizontal bands at M'^(rj7r) — 1 GeV^, there is a distinct 
enhancement at the highest three momenta in the bottom row of the figure. It 
turns out that this enhancement cannot be fitted adequately by constructive 
interference alone between two ao(980) bands; such an enhancement is not 
there at lower momenta. The enhancement in fact arises from a rapidly growing 
signal due to /o(2100), whose kinematic threshold is at 1525 MeV/c. This 
feature is important in two respects. Firstly it shows the presence of /o(2100) 
in rjr]. This was reported in Ref. [38]. Secondly the /o(2100) is so close to 
threshold that it appears only in the /o(2100)7r S-wave, which has = 0~ in 
pp. This observation helps identify the highest 0~ resonance. 

Table 4 shows the conspicuous resonant contributions in 777777°. The final col- 
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Fig. 37. Projections of r/r/vr data on to 7777 mass; histograms show the fit and beam 
momenta are indicated in MeV/c. 

umn shows changes in log Ukehhood when each resonance is omitted from the 
fit. All are statistically at least 14 standard deviation signals, though uncer- 
tainties in smaller amplitudes introduce systematic errors which reduce the 
significance level to 8-12 standard deviations. 

Singlet amplitudes are determined simply via Lcgcndre polynomials for pro- 
duction amplitudes. Angular distributions shown in Ref. [5] verify directly the 
presence of 0~, 2~ and 4~ states. There is evidence also in 37r° data for both 
0~ resonances, though they are determined best by the 777777 data. Both should 
be regarded as 2*. 

The 2~ state at ~ 1990 MeV is conspicious also in Stt^ data. Fig. 38(b) below. 
The combined fit to both channels gives an optimum mass of 2005 ±15 MeV 
and a width of 200±40 MeV, as shown below in Table 5. The observation of this 
2~ state in ao(980)77 is rather important. This resonance overlaps significantly 
with 712(1880) discussed below. However, 7r2(1880) has a very weak decay to 
Co (980)77 if ^11- The second order Legendre polynomial associated with 

the 2~ state at 1990 in 777777° data is very distinctive at 600, 900 and 1050 
MeV/c and leaves no doubt of its presence. This state should be regarded as 
3*. The 4~ signal is entirely consistent with a rather strong and distinctive 4~ 
signal in /27r in 37r° data and should be regarded as 3*. 

The 3"^ state at 2031 MeV is a confirmation of a very strong /27r signal in 37r° 
data, making this state 3*. The upper 3+ state is not well identified in 37r° 
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Name 




M 

(MeV) 


r 

(MeV) 


Decay 
Channel 


AS' 


vr 


0" 


2070 ± 35 


j-i on 

310t^°° 


ao(980)?? 
/o(1500)7r 


272 

324 


TT 


0" 


2360 ± 25 


3001100 


ao(980)r/ 
/o(2100)7r 


197 
434 




2- 


1990 ± 30 


200 ± 40 


ao(980)r/ 


531 


7r4 


4- 


2255 ± 30 


185 ± 60 


ao(980)r/ 


213 




2+ 


2255 ± 20 


230 ± 15 


/2(1980)7r 


363 


as 


3+ 


(2031) 


(150) 


02(1320)7? 
/o(1500)7r 


185 
164 


as 


3+ 


2260 ± 50 


2501^00 


oo(980)r? 


141 



Table 4 

Resonances parameters fitted to rj-qir^ data; values in parentheses are fixed from 
SyfO data. The final column shows changes in log likelihood when each channel is 
dropped from the fit and all others are re-optimised. 

data, so this should be regarded as 2*. The 2+ state at 2255 MeV confirms a 
state in Stt" and t^tt" and should be regarded as 3*. Its distinctive decay 
into /2(1980) will be discussed in detail in Section 10.4. 

6.3 Pp^Stt^ 

Table 5 shows resonance parameters from the better of the two solutions men- 
tioned above. This solution has log likelihood better by 1182. This is com- 
parable to the effect of removing one of the weaker resonances. The essential 
change in going to the poorer solution, shown in Table 6, is that the masses 
of the upper 4+ state and both of the upper two 2+ states slide downwards 
by ~ 40 MeV. There is also a small increase in the mass of the lower 4+ state 
from 2005 to 2030 MeV. What is happening is that the two 4+ states, which 
have similar values of r4, are tending to merge. Elsewhere it is observed that 
such merging of resonances of the same almost always gives small improve- 
ments in log hkelihood through interference effects. Generally this should be 
regarded with suspicion. The masses of the highest 4+ and 2+ states are 60-70 
MeV lower than for / = 0. If one takes the / = solution as a guide, this 
solution looks physically less likely, as well as having the poorer log likelihood. 

Despite differences of detail between the two solutions, the general pattern 
of masses and widths for 7 = and 7 = 1 is similar. It is noticeable that 
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* 


M 
(MeV) 


r 

(MeV) 


r 


AS 

(MeV) 


M{I = 0) 


Observed 
channels 


7r4 


4- 


3* 


2250 ± 15 


215 ±25 




11108 


2328 ± 38 




a4 


4+ 


3* 


2255 ± 40 


33o;^j° 


0.87 ±0.27 


2455 


2283 ± 17 


7rr/,7rr7',/27r 


04 


4+ 


2* 


20051^^ 


180 ± 30 


0.0 ±0.2 


2447 


2018 ± 6 


7rr7,/27r 


as 


3+ 


2* 


2275 ± 35 


350l^go 




3154 


2303 ± 15 


aori,f2V 


as 


3+ 


3* 


2031 ± 12 


150 ± 18 




18410 


2048 ± 8 


/o(1500)7r,a277,/27r 


a2 


2+ 


3* 


2255 ± 20 


230 ± 15 


-2.13 ±0.20 


2289 


2293 ± 13 


T^VJ2V, /2(1980)7r 


02 


2+ 


2* 


2175 ± 40 


310^90 


-0.05 ±0.31 


1059 


2240 ± 15 


/27r 


02 


2+ 


2* 


2030 ± 20 


205 ± 30 


2.65 ±0.56 


1308 


2001 ± 10 


vrr/,7r?7',/27r 


7r2 


2" 


2* 


2245 ± 60 


320tl°° 




2298 


2267 ± 14 


/27r,(cT?7) 


7r2 


2" 


3* 


2005 lb 15 


200 ± 40 




1633 


2030 ± 15 


ao??,/27r 


ai 


1+ 


2* 


2270^1^ 


305i^° 




2571 


2310 ± 60 


/27r 


77 


0- 


2* 


2360 ± 25 


300l^g° 




1955 


2285 ± 20 


ao??,/o(2100)7r,/27r 


TT 


0- 


2* 


2070 ± 35 


3101^°° 




1656 


2010^^^ 


ao??,/o(1770)7r,/27r 


ao 


0+ 


2* 


1995 ± 30 


300 ± 35 




2134 


2040 ± 38 


7rr7,7rr7',77(1440)7r 


02 


2+ 




(1950) 


(180) 


-0.05 ± 0.30 


2638 


1934 ± 20 




vr2 


2- 


4* 


(1880) 


(255) 




7315 


1860 ± 15 




ai 


1+ 




(1930) 


(155) 




2609 


1971 ± 15 




01 


1+ 




(1640) 


(300) 




232 






vr 


0" 




(1801) 


(210) 




2402 







Tabic 5 



Masses and widths of fitted I = 1, C = ±1 resonances. Values in parentheses are 
fixed from other data. Entries in the lower half of the table are below the avail- 
able mass range. Column 7 shows changes AS in log likelihood for 37r° when each 
resonance is removed from the fit and all remaining parameters are re-optimised. 
Column 8 shows masses for I = 0, C = +1 resonances from Table 2 for comparison. 

most 1 = 1 masses tend to lie 20 MeV lower than for / = 0. Since the poorer 
solution tends to drag the masses of the upper 2+ and 4+ states down, there is 
the possibility that this ambiguity is everywhere having the effect of lowering 
masses, through correlations between partial waves. 

One of the problems in the analysis is that angular distributions of 3-body data 
change rapidly between the two lowest beam momenta, 600 and 900 MeV/c. 
In order to accomodate this feature, both the lower 3+ and 4+ states must be 
rather narrow. Popov [57] reports an 7 = 1, = 3+ peak in E852 data for 
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Name 




M 


r 


AS 


r 






(MeV) 


(MeV) 






04 


4+ 


2220 ± 20 


345 ± 65 


2852 


0.30 ±0.31 


04 


4+ 


2035 ± 20 


135 ± 45 


2063 


0.0 ±0.2 


02 


2+ 


2235 ± 35 


200 ± 25 


3658 


-1.74 ±0.36 


a2 


2+ 


2135 ± 45 


3051^0 


1685 


-0.56 ± 0.53 



Table 6 

Masses and width of resonances in the alternative solution. The last column shows 
changes AS in log likelihood when each resonance is removed from the fit and all 
remaining parameters are re-optimised. 

[6i(1235)7r] at ~ 1950 MeV, but fits to a much lower mass of M = 1873±1±18 
McV with r = 259 ± 18 ± 21 McV. He also reports an 04 at 1952 ± 6 ± 14 
MeV in [u;p]l=2 with T = 231 ± 11 ± 21 MeV. 

Fig. 38 shows as full curves the intensities fitted to individual /27r partial 
waves in the better solution. Intensities for single energy solutions are shown 
by black squares. This gives some idea of the stabihty of the fit. 

Because of the ambiguity between two solutions, none of these states can be 
regarded as 4*. Both 4+ states are certainly required by nr) and /27r data; the 
upper one is required also by iri]' data. The lower 3"'" state is exceptionally 
strong, Fig. 38(g) and must be regarded as well established. The 2+ states are 
poorly separated between and ^^2, but a solution is not possible without 
three of them in the mass range 2000 MeV upwards. Except for the /2(2255), 
which appears distinctively in 777777, they should be regarded as 2* because of 
the ambiguity in the solution. 

The upper 2~ state is required but poorly determined and can only be classed 
as 2*. A lower 2~ state is definitely required by the peak in /27r in Fig. 38 (b). 
Recently the E852 collaboration has presented evidence for a similar resonance 
in /i(1285)7r with M = 2003 ± 88 ± 148 MeV, F = 306 ± 132 ± 121 MeV [58]. 

For = 1+, both /27r data and 777777 require a 1+ state around 2270 MeV, but 
it can only be regarded as 2*. There is a strong 1+ intensity at 600 MeV/c. 
The 0~ states are better identified by 777777 data. 

As remarked above, there is a large change at the lowest momenta between 
the earhest analysis of /27r [53] and the final analysis of Ref . [4] . The problem 
is an ambiguity between ^Pi, and ^^4 amplitudes and the rapid change of 
angular distributions with mass. The ambiguity between ^Pi and ^P2 is the 
principal problem, so both can only be regarded as tentative identifications 
at present: something to watch for in other experiments. The very strong 
state is well estabhshed independently of the others. Good information from 
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Fig. 38. Intensities of individual /27r partial waves for the preferred solution. Squares 
and triangles show free fits at single energies. Full curves are for /27r final states. 
Dotted and chain curves in (a) show respectively the 4~ /27r amplitude with L = 5 
(triangles) and the gtx contribution. Dotted curves in (c) and (p) show o-tt contri- 
butions. In (p), points and the full curve show integrated cross sections for Zti^ . 

any other source resolving the ambiguity w^ould greatly stabilise the complete 
solution for / = 1, C = +1. One should be alert for a 2+ state decaying 
dominantly to pcu, as partner to /2(1910 — 1934) — > lolo. The strong rise in 
Fig. 38(i) at the lowest momentum suggests a contribution from a low mass 
state around 1950 MeV. However, this is right at the bottom of the available 
mass range and cannot be regarded as established at present. 



6.Ji. Evidence for 772(1 880) 



Evidence was presented for 7^2(1870) in Section 4.3. An I — 1 partner is to 
be expected and is indeed observed in pp — > 7r27r, 7r2 — > 02(1320)77, i.e. in 
the TjriTi^Ti^ final state [59]. The entire pattern of production cross sections is 
remarkably similar to that for rjir^Tr^ir^. 
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Fig. 39. Full histograms show fits to rjrjTr^Tr^ data at 1200 MeV/c compared with 
data; dotted histograms show fits without 7r2(1880); (a) the rjrj-K mass distribution, 
(b) the TTT] mass distribution, (c) the rj-qn mass distribution for events with a ttt] 
combination within ±65 MeV of 02 (1320), (d) and (e) nr] and r]r] mass distributions 
for events with M{r]r]'K) 1850-2050 MeV. Units of mass are MeV. 

The evidence for 7r2(1880) is found at beam momenta 900, 1050 and 1200 
MeV/c (as for r/2(1860)). Fig. 39 shows data at 1200 MeV/c. In (a), dashed 
histograms show the fit without 7r2(1880). There is a hint of a peak in ryr^Tr. 
In (c), there is clear evidence that this is associated with the 02(1320). The 
essential evidence comes from (c), (d) and (c). Without 7r2(1880), the dashed 
histogram in (c) fails to fit the observed 02 (1320) signal. Note that the normal- 
isation of the curve cannot be scaled up arbitrarily to fit the data, since it is 
governed by fitted amplitudes. In (d), the 02 signal in the T^r^Tr mass range 1850 
MeV (approximate threshold) to 2050 MeV is clear. In (e) the rjrj enhancement 
below 1350 MeV arises as a reflection of the decay channel 02(1320)77. 

Optimum parameters are M = 1880 ± 20 MeV, F = 255 ± 45 MeV. The 
signal cannot be explained as the high mass tail of 772(1660): it is too strong 
and the 7r2(1670) is too narrow to fit observed cross section. From Fig. 39(a), 
the signal clearly does not arise from r^r^vr masses near 2005 MeV. The angular 
distribution for decay to 02(1320)7; is consistent with S-wave decay, as expected 
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Fig. 40. Data of Daum et al. [62]. The upper figure shows the /2vr D-wave intensity; 
the lower figure shows the phase difference between 7r2(1670) and [/27r]L=2. 

so close to threshold; it is inconsistent with L — 1 decays. 

The VES collaboration also reports a threshold enhancement in a2(1320)r; with 
— 2~ in their rjrjn' data [60] and also in jyTr+Tr^Tr^ data where 02(1320) — > 
7r+7r~7r°. They observe a strong peak at the same mass in [/2(1270)7r]L=2 [61]. 
Daum et al. [62] also observed a peak at 1850 McV in [/2(1270)7r]i,=2 with a 
width of ~ 240 MeV. Their result is reproduced in Fig. 40. There is a clear 
peak and a phase motion with respect to 7r2(1670) entirely consistent with a 
second 2~ resonance associated with the peak. However, they interpreted the 
data conservatively as arising from interference between 772(1 670) and a higher 
772(2100); that interpretation fails for Crystal Barrel data. 

Popov [57] reports a 2^ peak in E852 data for [ci;p]l=i for the combination of 
u and p spins = 1. He quotes M = 1890 ± 10 ± 26 MeV, F = 350 ± 22 ± 55 
McV. The mass agrees well with Crystal Barrel; the width is slightly higher. 
Eugenio [63] reports a similar 2^ peak in E852 data for the ririn^ final state. In 
view of all these observations, the 772(1880) should be regarded as a 4* state. 

In section 5.3, it was reported that 772(1870) has a branching ratio to /o(1270)?7 
much stronger than to 02(1320)77 and inconsistent with SU(3) for a qq state. 
The 772(1880) is observed strongly in decays to 02(1320)77 by VES but no decay 
to the /277 S-wave is observed; although no actual branching ratio is reported 
at present, it seems likely that 772(1880) is also inconsistent with a qq state 
and SU(3) relations for decays. 



7 / = 1,C=— 1 resonances 

Data are available from 77" 77+, un'^ and curjir'^ final states. The latter two 
channels have been measured for uj — > 77''7 at all momenta and for ou — > 
77+77-77° at 6 momenta: 600, 900, 1200, 1525, 1642 and 1940 MeV/c. 
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The analysis of pp — > 7r~7r+ is reported in Ref. [13]. The I — 1 states with 
J^'-^ = 1 , 3 and 5 cause forward-backward asymmetries in both da/dQ 
and Pda/dQ via interferences with / = states; these asymmetries are con- 
spicuous and lead to good determinations of C = — 1 amplitudes. A particular 
feature at low momenta is a very strong p3(1982) which interferes with /4(2050) 
over the whole momentum range around 750 MeV/c. Its mass and width are 
determined very accurately by the 7r~7r"'" data. This resonance then serves as 
a useful interferometer in ujti^ and uriiT^ data in the low mass range. A second 
J^*-" = 3 state is required also at 2260 MeV and serves as a similar inter- 
ferometer at high masses. The p5(2295) appears as a large signal in tt'tt'^, see 
Fig. 14 of Section 5.1. 

Attention needs to be drawn to one point concerning the 1 state at 2000 
MeV. The best mass determination of this state comes from 7r~7r+, because 
of the PS172 data available in small momentum steps down to a mass of 1912 
MeV. As emphasised above, the combination of differential cross sections and 
data from a polarised target determines phases well for all partial waves. The 
Argand diagram for ^Di on Fig. 23 of Section 5.4 requires a phase variation 
approaching 180° between 1900 and 2100 MeV, hence a 1 resonance. The 
^Si amplitude shown on the same figure has a 1/v behaviour close to the pp 
threshold. When one factors this 1/v dependence out, the ^5*1 amplitude also 
requires the same 1 resonance. For ujtt^ and ujrin^ data discussed below, 
data from a polarised target are not available, so determinations of its mass 
and width are less reliable. Once it is established, the 1 amplitude again 
serves as a useful interferometer in utt^ and curjTr^ data at low masses. 

Analysis of pp un^, u — > n^j is reported in Ref. [64]; ujrjTr^ data with 
uj — > 7r°7 are reported in Ref. [65]. Although this uj decay mode is identified 
cleanly, polarisation information for the cu is carried away by the photon, whose 
polarisation is not measured. This was the reason for studying cu Tr+Tr"?!^ 
[66]; the combined analysis of all / = 1, C = —1 channels is reported there. 
The matrix element for u Tr+Tr^Tr" is relativistically ^^/s-ysP^P^P^ where p are 
4-momenta of decay pions. Non-relativistically, this matrix element becomes 
Pi A pj in the rest frame of the cu, where pi are 3-momenta of any two pions. 
The polarisation vector of the u is therefore described by the normal n to its 
decay plane. A measurement of this decay plane provides complete information 
on the polarisation of the uj. That improves the amplitude analysis greatly. 
Conclusions should therefore be taken from Ref. [66]. 

The measurement of both decay channels also provides a valuable experimental 
cross-check, illustrated in Fig. 41. Angular distributions are measured mostly 
in the backward hemisphere for the u, because of the Lorentz boost from 
centre of mass to the lab. Figs. 41 (e)-(j) show as points with errors the angular 
distributions for a;7r° from cu — > 7r'*'7r~7r° decays. The corridors show angular 
distributions from the partial wave analysis of a; ^ 7r°7. Both are corrected 



56 




800 



700 800 900 




300 - (e) 600 MeVA 



CO 

c 

> 200 fi 



100 - 



-1 -0,5 0.5 



800 



0.5 0.5 



400 - 



200 



-0.5 





MeV 


300 


400 


; (c) J^- 


200 




000 




200 


800 




1 50 


600 
400 




100 


200 




50 











700 800 900 

MeV 



0.2 0.4 0.6 0.8 1 




1000 




(i) 1642 MeV/c 




0.5 1 
cos 




Fig. 41. M(7r+7r^7r°) from the data selection for (a) ujir^ data, (b) iorjir^ at 900 
MeV/c; the number of events v. the square of the matrix element for u decay in 
(c) a;7r°, (d) wqir^; (6)-(j): curves show the error corridor for uir^ differential cross 
sections where u decays to they are corrected for angular acceptance; points 
with errors show results for u) decays to Tr+vr^vr^; the vertical dashed lines show the 
cut-off which has been used. 
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for acceptance. The observed agreement is a major cross-check on cahbrations 
of the detector, data selection and Monte Carlo simulations. All of these are 
completely different between all-neutral data and charged data. One aspect of 
this is the treatment of the vertex. For charged data, the vertex is determined 
by the intersection of the two charged tracks. For neutral data, the vertex is 
instead taken to be at the centre of the target. This assumption distorts the 
kinematic fit slightly, but the correction is only 1% to da/dfl for all-neutral 
data. This procedure is safer than fitting the vertex freely; that alternative 
leads to strong and dangerous variations of acceptance with cosO^^. 

Backgrounds are clearly visible in Figs. 41 (a) and (b) under the uj. The 
backgrounds are checked in (c) and (d). There, the number of events is plotted 
against the square of the matrix element for u decay. The linear dependence 
checks that the u; signals are clean. The intercepts at |Mp = check the 
magnitudes of the backgrounds. 



7. 1 Vector polarisation of the oo and its implications 



A remarkable feature is observed for u polarisation, with major physics im- 
plications. Both vector polarisation Py and tensor polarisations T20, T21 and 






Fig. 42. Vector polarisation Py for u;7r° data at (a) 900, (b) 1200, (c) 1525 and (d) 
1940 MeV/c compared with the partial wave fit (histogram); (e)-(h) Re T21 at the 
same momenta. 
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Fig. 43. (a)-(d) Re T22 and (e)-(h) T20 for ujtt^ data at 900, 1200, 1525 and 1940 
MeV/c, compared with the partial wave fit (histograms). 

T22 are measured. For details and formulae, see Ref. [66]. Vector polarisations 
are shown in the upper half of Fig. 42. They are consistent with zero almost 
everywhere. Large tensor polarisations are observed at all momenta in the 
lower half of Fig. 42 {Re T21) and in Fig. 43 (T20 and T22). Vector polarisation 
depends on the imaginary parts of inteferences between partial waves. Ten- 
sor polarisations depend on moduli squared of amplitudes and real parts of 
interferences. In particular. Re T21 contains real parts of precisely the same 
interference terms as appear in Py- Formulae are given by Foroughi [67] for 
pp — > dir'^, but the same formulae apply to pp — > a;7r°. The implication of zero 
vector polarisation is that relative phases of amplitudes are close to or 180°; 
the physics implications will now be discussed. 

Consider first a given such as 1+. This partial wave decays to con with 
orbital angular momenta L = or 2. It is not surprising that partial waves 
to these final states have the same phase, since a resonance implies multiple 
scattering through all coupled channels. Indeed, it is found that relative phases 
are consistent with zero within errors for all resonances. [The final analysis 
can then be stiffened by setting the relative phase angle to zero.] 

Vector polarisation Py can arise from interference between singlet waves 1"*", 
and S"*". Likewise, it can arise from interference between triplet waves 1~, 
2~, 3~, 4~ and 5~. It is remarkable that these further interferences also lead 
to nearly zero polarisation. It implies coherence between different J^. There 
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is no full explanation for this result, but an interpretation is suggested here. 



The incident plane wave may be expanded in the usual way in terms of Leg- 
endre functions: 



gifc. ^^^2£+l)i^P^(cos^). 



(20) 



As a result, all partial waves of the initial state are related in phase through 
the factor i^. Fig. 44 shows Argand diagrams for all partial waves. Let us take 
as reference for triplet states the large ^G4{£ = 3) amphtude. If one compares 
this by eye with diagrams for and ^D^, it is clear that there is a phase 
rotation of the diagam by ~ 180° for "^^2; is similar in phase to ''G'4. The 
Argand loop for ^5*1 is again rotated by ~ 180° with respect to '^(^4. Singlet 
and triplet states do not interfere in da/dQ, Py, T20, T21 and T22. Therefore 
the phase of singlet states with respect to triplet are arbitrary in Fig. 44. The 
Argand loops for ^Pi and ^F^ are broadly similar, though there is a large offset 



< 




Re A 

Fig. 44. Argand diagrams for partial waves fitted to llitt^ in the combined analysis 
with corjTT^ and 7r"'"7r~. Crosses show beam momenta 600, 900, 1050, 1200, 1350, 
1525, 1642, 1800 and 1940 MeV/c; all move anti-clockwise with increasing beam 
momentum. L is the orbital angular momentum in the ujtt channel. 
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Fig. 45. Projection of lot]?:^ data on to M^-K^rj); full histograms show the fit. 
in the ^Pi L — amplitude. 

A similar phase coherence is reported for pp — > 7r~7r+ [1]. There, both 7 = 
and I — 1 amplitudes are present. The amplitudes which differ by 1 in orbital 
angular momentum for pp are ~ 90° out of phase, as suggested by eqn. (21). 

An analogy is given in Ref . [66] . Suppose a bottle is filled with irregular shaped 
objects, such as screws. If the bottle is shaken, the objects inside settle to a 
more compact form. Returning to resonances, they have well defined masses, 
widths and coupling constants which make them individual. But a suggestion 
is that, in the interaction, phases adjust to retain as much coherence as possible 
with the incident plane wave. This is a natural way of inducing order between 
the observed resonances. It helps explain why all Argand diagrams show a 
smooth variation of phases with mass. 
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M(pp) (GeV) 



Fig. 46. Integrated cross section for pp — > wTyvr compared with the fit. It is corrected 
for acceptance and for branching ratios of a; — > tt^j, r/ — 77 and 7r° — 77. 

7.2 007)77^ data 

Fig. 45 shows the tt^t) mass projection from curjir'^ data. There are clear peaks 
due to ao(980) and 02(1320). The strong increase in the 02 signal at the higher 

masses is due to the opening of the a2UJ threshold at 2100 MeV. This threshold 
is also observed through its direct effect on the integrated cross section, shown 
in Fig. 46. In contrast, integrated cross sections for most other channels peak 
at low mass, e.g. Fig. 10(a) for Tr-Tr"^ in Section 3. 

A feature of the urin^ data is shown in Fig. 47. There is a clear peak at 
1650 MeV in ujt). The amplitude analysis favours J^*" = 1 strongly over 
jPC ^ ^+-^ 2— or 3—. The fitted mass and width M = 1664 ± 17 MeV, 
r = 260±40 MeV agree well with PDG averages for u;i(1650). The final state 
a;i(1650)7r contributes strongly to the l+~ (^Pi pp) resonance at 1960 MeV. 

7.3 Observed resonances 

Table 7 shows resonances from the final combined analysis of all channels. The 
complete spectrum of expected qq states is observed except for the ^Gs state 
and the highest ^5"! state expected near 2400 MeV. The state may easily 
be hidden close to the strong upper state; indeed, that resonance has a 
ratio r3 of ^G^/^D^ amplitudes significantly different from 0; this hints at the 
presence of either a ^G^ state or quantum mechanical mixing with it. There 
is even a hint of the last ^Si state. A scan of the mass of the 1~ state at 2265 
MeV reveals two distinct optima in log hkelihood, one at 2265 MeV and a 
second weaker one at 2400 MeV; however there is not sufficient evidence to 
isolate this state so close to the top of the available mass range. 
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600 MeV/c 



900 MeV/c 




Fig. 47. Projection of uotj-k^ data on to M{r]io); full histograms show the fit including 
an a;(1665) with T = 260 MeV. 

Intensities of partial waves fitted to a;7r° are shown in Fig. 48. A remarkable 
feature is a very strong state at 2230 MeV. In all other data, G-states 
are weak, presumably because of the strong £ — A centrifugal barrier in the pp 
channel. The analysis has been checked in many ways to ensure that there is 
no cross-talk with ^D2, which is quite weak at high momenta. As an example, 
analyses have been carried out at single energies. The fitted ^G^ intensity 
follows the behaviour expected from the centrigugal barrier. It falls close to 
zero at low momenta where the intensity is large; this is illustated in Fig. 
4 of Ref. [64]. So there is no doubt that the huge ■'(^4 amplitude is a real effect. 
There is no obvious explanation. It must indicate that for some reason this 
state decays dominantly through ujtt. The low mass of the state is also 
surprising; it is lower than the nearby F-states. This is in contrast to the well 
estabhshed behaviour of / = states where masses of G-state resonances are 
higher than those of F-states. 

Table 7 shows in the right-hand entries the firmly established decay channels 
where each resonance is observed. The star rating depends on the number of 
channels in which each resonance is observed and the stability of resonance 
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Table 7 

Resonance parameters from a combined fit to ujtt^, urjii^ and 7r~7r'^, using both 
u> TT^j and uj vr+vr^vr'^ decays. Values of r are ratios of coupling constants 
gj+i/gj-i for coupling to pp or utt^. Columns 6 and 7 show changes in = log 
likelihood when each resonance is omitted and others are re-optimised. In the final 
column, uj' denotes a;i(1650); channels in parentheses are weak. 

parameters. Figs. 49 and 50 show triplet and singlet partial waves for pp — > 
curjTT^ from the final analysis of Ref. [66]; these figures are an update of those 
from Ref. [65], but changes are small. Brief comments will now be given, 
starting with the highest spins. 

The ^5(2300) is already known from GAMS data [68]. The determination from 
pp TT~TT^ is more accurate, from the large 5 amplitude observed there. 
GAMS omitted any centrifugal barrier for the production process; this barrier 
lead to a peak position in Crystal Barrel data of 2340 MeV in agreement with 
the mass of 2330 ± 35 MeV quoted by GAMS. 

The 4r state at 2230 MeV is observed independently in ujii decays with L = 3 
and 5 and also in ao(980)a; in Fig. 49. Since it is the dominant signal at high 
masses in lot:, it is very clearly established and classified as 4*. 

The upper 3~ state is observed in 5 channels; the lower state at 1982 MeV is 
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Fig. 48. Intensities of partial waves fitted to lott data; L is the orbital angular 
momentum in the decay. 

strong in 7r^7r+ and accurately determined. Both are therefore classified as 4*. 
Indeed, the trajectory of 3 states {pp ^D-^) in Fig. 3(c), including the very 
well known p3(1690), provides one of the best determinations of trajectory 
slopes. 

The upper 2~ state at 2225 MeV is not given 4* status because of correlations 
with 3^ and 1^ partial waves. The lower 2^ state at 1940 MeV is one of the two 
most weakly established states. However, its Argand plot on Fig. 44 requires 
a rapid phase variation at low mass because the amplitude must go to zero at 
the pp threshold. So a resonance is required, even if it is at the bottom of the 
available mass range and therefore somewhat uncertain in mass. 

All three 1 states are observed in three channels. The p(2150) was found 
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Fig. 49. Intensities of dominant triplet partial waves in ujrjTr data; dashed lines refer 
to the second title in panels. 

by GAMS [69,70] and now deserves 3* status. The mass scan for the upper 
one gives a distinct optimum at 2265, but some indication of a second at 2400 
MeV. It is possible that there are two states in the high mass range and they 
cannot be fully resolved. The state at 2265 MeV is therefore classified as 2*. 

Now consider singlet states. The 3"' state C'Fs) at 2245 MeV is observed in 
4 channels; intensities of singlet partial waves for uriTi^ arc shown in Fig. 50. 
Its presence is required in utt by the 360° loop in the Argand diagram of Fig. 
44. However, there is a rather large error of ±50 MeV in its mass. The lower 
3+~ state is the strongest contribution to lot: and deserves 4* status. 

The partial wave analysis requires a rather large 1"*"" intensity in the a;(1650)7r 
S-wave, as shown in the bottom left panel of Fig. 50. It peaks above 2 GeV 
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Fig. 50. Intensities of singlet partial waves in cjrjTr data; dashed lines refer to the 
second title. 

because of the i = 1 centrifugal barrier for production and also the increasing 
phase space for production of C(j(1650). It is identified directly by the peak 
in the urrj mass distribution. The ujtt^ data also require a rather strong ^Pi 
intensity at low masses. It cannot be eliminated without a large change of 289 
in log likehhood. The Argand loop of Fig. 44 for lott^ requires a strong phase 
variation when one remembs that the amplitude must go to zero at the pp 
threshold; this demands a low mass l"*" resonances, even though there is some 
uncertainty in its precise mass. 

A final remark is that Matveev reports a new b'l ed, M — 1620 ± 15 MeV 
with r = 280 ± 50 MeV at the recent Aschaffenburg conference [71]. This is 
included on Fig. 3(a). 
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8 I — 0, C — —1 resonances 



There are three pubhcations of Crystal Barrel data with these quantum num- 
bers. The first was for a;7r°7r° data with u; decaying to 7r°7 [72]. This was 
followed by an analysis of ujr] where uj again decays to 7r°7 [73]. The final 
combined analysis added data for both channels where u tt^tt^tt^ [74]. 
The measurement of uj polarisation in these data improves the stability of the 
analysis considerably, so conclusions will be taken from that work. 

Features of the data are very similar to those for / = 1, C = — 1 and need 
little comment. Figs. 51(c)-(j) compare angular distributions for uji] in decays 
via n^j (shown by the corridors) and via u tt+tt^tt^ (points with errors). 
There is agreement within errors, as there was for a;7r°. This checks the data 
selection and Monte Carlo simulations, but with lower accuracy than for ujii'^ 
because of lower statistics. Fig. 52 shows vector and tensor polarisation Re T21. 
Other tensor polarisations are large, as for ujti^ , and are shown in Ref. [74]. 
Again the vector polarisation is consistent with zero, indicating the same phase 
coherence amongst the amplitudes as for ujt:. 
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Fig. 51. M(7r^7r^7r'^) from the data selection for (a) corj data, (b) lott^tt^ at 900 
MeV/c; the number of events v. the square of the matrix element for u decay in 
(c) LOT], (d) WTT^TT^; (e)-(j): curves show the error corridor for urj differential cross 
sections where uj decays to they are corrected for angular acceptance; points 
with errors show results for co decays to Tr+vr^vr^; the vertical dashed lines show the 
cut-off which has been used. 
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Fig. 52. Vector polarisation Py for ut] data at (a) 900, (b) 1200, (c) 1525 and (d) 
1940 MeV/c compared with the partial wave fit (histogram); (e)-(h) Re T21 at the 
same momenta. 

Fig. 53 shows intensities of partial waves fitted to oorj and Fig. 53 shows Argand 
diagrams. Fig. 55 shows intensities of the large partial waves fitted to lutt^tt^ 
data. 
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Fig. 53. Intensities of partial waves fitted to iorj data; L is the orbital angular 
momentum in the decay. 

The analysis of this family of / = 0, C = — 1 states is the weakest of all four 
families. Statistics for the u>r) channel are only ~ 5000 events per momentum. 
In the Lon^n^ data, there is a strong /2(1270) signal but it is superimposed 
on an cua physics background. Uncertainties in parametrising the broad tttt 
S-wave amplitude above 1 GeV introduce possible systematic errors into this 
physics background and, via interferences, into the determination of resonance 
masses in the uf2 channel. However, without the ua channel, the data cannot 
be fitted adequately, so there is no choice but to include it. 

There is evidence in the final analysis for all expected qq states except for two 
1^. However, several of them have a 2* rating because of (i) low statistics for 
CUT}, (ii) the difficulties caused by the ura physics background in a;7r°7r°, (iii) the 
lack of an accepted interferometer to establish phase variations. Table 8 shows 
fitted resonances and their star ratings. None can be classified as 4* because 
of the limited number of decay channels. 

It is best to start with the low mass range and triplet states. The ^D^ partial 

wave in lvtt'^'k^ data of Fig. 55(d) is very strong in the final state down to 
the lowest beam momentum. This cannot be confused with a threshold effect 
in ^Si because ^5*1 and '^^3 give rise to quite different tensor polarisations 
of the /2, hence distinctive decay angular distributions. For cu — > 7r''"7r~7r° 
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Fig. 54. Argand diagrams for partial waves fitted to ujrj in the combined analysis 
with u-K^-K^. Crosses show beam momenta 600, 900, 1050, 1200, 1350, 1525, 1642, 
1800 and 1940 MeV/c; L is the orbital angular momentum in the wq channel. 

decays, the correlation between the normal to the uj decay plane and the 
decay angular distribution of the /2 is very sensitive to the different Clebsch- 
Gordan coefficients for ^-Da, '^-02, '^Di and amplitudes. The conclusion 
is that the selection of the ^D^, initial state is very distinctive. This can be 
demonstrated by running a crude initial analysis at low momenta with just one 
partial wave at a time. A small but significant 6i7r D-wave amplitude provides 
further evidence for the ^D^ initial state. 

Results for ujt] are similar. Again the normal to the uj decay plane in a; — > 
7r"'"7r~7r° decays provides good separation between ^D^, and ^Di. The 
momentum dependence separates pp S and D-waves. The differential cross 
sections for uji] look like (1 + 2cos^^) after correction for acceptanc; here 6 
is the production angle. The strong cos^ 6 term requires initial S or D-states 
(and final P-states). The data cannot be reproduced by singlet partial waves. 
The ^Pi amphtude gives the isotropic term. In the earliest analysis of Ref. 
[73], the amphtude was neghgible. 
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Fig. 55. Intensities of the large partial waves fitted to uiir it data. 

Argand diagrams of Fig. 54 for ^D^, and "^Di lot] final states require phase 
variations approaching 180° between the pp threshold and 2100 MeV (a beam 
momentum of 1050 MeV/c). These partial waves dominate the tor] differential 
cross section up to 1200 MeV/c. The rapid decrease with momentum of the 
ujr] cross section forces all of the partial waves to agree with quite narrow 
widths: 115 MeV for ^Dg, 175 McV for and 195 MeV for ^Di. It is not 
possible to make them wider and fit the observed cross sections. Observed 
phase relations between ^^3, and ^Di make resonances in all three partial 
waves unavoidable in the same mass range. The ^5*1 partial wave is difficult 
to separate from a pp threshold effect, but its phase falls naturally into line 
with that observed for ^D^. 

In view of the well established I —1 ^D^ state at 1982 MeV, these results are 
not surprising. 

An isotropic component of the wq differential cross section is fitted naturally 
by ^Pi; again the normal to the uj decay plane in a; — > 7r'*"7r~7r'^ isolates this 
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Table 8 



Resonance parameters from a combined fit to torj and ton n using both lo tt 7 
and uj TT^n^TT^ decays. Values in parentheses are fixed. Values of r are ratios 
of coupling constants gj^i/gj^i for coupling to pp or cor]. Columns 6 and 7 show 
changes in S = log likelihood when each resonance is omitted and others are re- 
optimised. 

partial wave. A difficulty here is that there is no interferometer to define the 
phase of the ^Pi amplitude. However, its mass dependence fits naturally to 
a Breit-Wigner amplitude multiplied by the factor ppp/p of eqn. (3), section 
4.1. The charged decay mode of the u; also successfully locates a low mass ^F^ 
state at 2025 MeV; this state evaded detection in the analysis of a; — > 7r°7. 

However, it is clear that all these resonances arc correlated in masses and 
widths. There is no established resonance to act as interferometer, so resonance 
masses and widths are derived from the simple assumption of eqn. (3) for Breit- 
Wigner resonances. Results for masses come out in fair agreement with 1=1, 
C = — 1; there the p3(1982) is derived precisely from interferences with the 
known /4(2018). In view of the assumptions in the / = 0, C = — 1 analysis, it 
is not surprising that / = S D-states have masses which are systematically 
slightly lower than 7 = 1, C = — 1. This may not be a real effect. 

The conclusion is that all of the strong 1~, 2~ and 3^ states at low masses 
make a pattern consistent with / = 1; they deserve 3* states for their existence, 
subject to the reservation that there may well be a systematic error in their 
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mass scale. The weak ^Pi and states at low mass are classified as 2*. 

Next consider high mass states. The Argand loops of Fig. 54 for ^Di, 
and ^D^ all move through approximately 360° between the pp threshold and 
the highest mass, 2410 MeV. This is consistent with a group of higher 1~, 
2~ and 3~ states. The structure in the Argand loop for ^D^ indicates two 
separate states. This structure is apparent also in the ^D^ intensity of Fig. 53. 
Intensities of ^Di and partial waves drop to low values at high masses in 
Fig. 53, but the loops in Argand diagrams are distinctive. 

The situation concerning the 3~ amplitudes at high mass is somewhat sur- 
prising. The measured tensor polarisation in the 00/2 final state in o^tt^tt^ data 
appears to require a strong peak in Fig. 55(g). The data fit naturally to a 
strong ^6*3 state at 2255 MeV and a weaker ^D^ state at 2285 MeV. However, 
this result should be treated with caution. The polarisation of the uj in the final 
state is sensitive to the orbital angular momentum L in the final state. How- 
ever, it is much less sensitive to the angular momentum £ in the production 
process; ^G^ and ^D^ states differ only through Clebsch-Gordan coefficients 
for couphng of pp helicity states. To be sure of the separation between ^^3 
and ^Gs, data from a polarised target are needed and are not available at 
present. A balanced view of the results is that there is evidence for both ^G^ 
and ^1^3 states, but they arc not well separated by present data. In Table 8, 
they are shown as 2*/3*; the 3* refers to their existence, but the 2* refers to 
the imprecise identification of masses. 

The upper 2 state is weak in ujrj and determined only through its interfer- 
ence with other neighbouring states and with the lower 2^ state. It is therefore 
only 2*. 

The upper ^Di state at 2205 MeV appears weakly in uf] and is determined 
mostly by cua data. Because of the uncertainties in parametrising the broad 
(7, it can be classified only as 1*. For / = 1, there are two neighbouring states 
at 2110 and 2265 MeV. This suggests that the 7 = state may be a blurred 
combination of two states. 

The 5~ state at 2250 MeV appears quite strongly in both ujt] and irbi, Fig. 
55(g), but the mass has a rather large error, so it only deserves 2* status. The 
4~ state appears only in cor] and is 1*. 

Turning to singlet states, the ^F^ resonance at 2275 MeV appears strongly in 
U!f2 and is also required by ur] with L — A. The Argand diagram for ^^3 with 
L — A completes a full 360° loop which can only be fitted with two 3+ states. 
Since the mass is well determined, this state is given 3* status. The upper 1"*" 
state appears weakly in cur] and is required by a;/2. It deserves 2* status. 

In summary, the pattern of/ = 0,C = — 1 states is very similar to that of 
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I — 1, but is less well defined because of (i) the lack of data from a polarised 
target, (ii) the uncertainties in fitting the ua channel. Several of the states 
lying at the bottom of the mass range may be subject to systematic errors in 
mass. That situation is better for / = 1 because of the strong 3~ state at 1982 
MeV. This state is accurately determined by tt^tt"*" data and then acts as an 
interferometer with other partial waves. No such interferometer is available 
for 7 = 0, C = -1. 



9 Comments on systematics of qq resonances 

The first obvious comment is that there is no proof that these are qq states 
rather then hybrids. The regularity of the spectra suggests an origin common 
with accepted qq states at low mass. There is firm evidence for just one hybrid 
with exotic quantum numbers J^^ — from the E852 collaboration at 1593 
MeV [75] . There is also evidence for structure with the same quantum numbers 
near 1400 MeV from the E852 collaboration and from Crystal Barrel data at 
rest [76,77]. Hybrids are expected closer to the 1800 MeV mass range. So it 
seems possible that the 1400 MeV structure is connected with the opening of 
the first S-wave inelastic thresholds with these quantum numbers: 6i(1235)7r 
and /i(1285)7r. Further hybrids are expected in the mass range above 1800 
MeV, but they are predicted to be broad and may be difficult to locate. 

9. 1 Slopes of trajectories 

Figs. 1-3 of Section 1 summarise the observed spectra. Straight lines are drawn 
through them. An obvious comment is that this may well be an approxima- 
tion. It is to be expected, for example, that the opening of strong thresholds 
may perturb the precise masses of resonances. An obvious example is /2(1565), 
which sits precisely at the ujuj threshold and couples strongly to that chan- 
nel [78]. Therefore Figs. 1-3 should be taken as a general indication of the 
systematics. The straight lines do have predictive power. For example, the 
6'^ (1620) was actually predicted very close to this mass from the known mass 
of 6i(1230) and was found just where expected [71]. A second similar example 
is the recently discovered /ii(1595). 

There are still some missing states. In the mass range 1600-1800 MeV, an /i 

is expected around 1670 MeV and an rj close to 1680 MeV. These low spin 
states are difficult to identify in production processes in the presence of strong 
signals of higher spins. The ^D^ states of both isospins are also missing around 
1660 MeV. The 0+ sector for 7 = 1 is poorly represented at present. Further 
states are expected around 1760 and 2380 MeV; the latter is right at the top 
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Table 9 

Slopes of trajectories. 

of the Crystal Barrel range. Also the whole situation concerning ^5*1 states is 
obscure at present. 

Table 9 summarises the slopes of/ = 0,C = +l resonances, which are the 
best established; it also includes a very well determined slope from the / = 1, 
J^^ — 3 family, which have accurately known masses. The mean slope is 
1.143 ±0.013 GeV^. At high masses, it is expected that spin-orbit, tensor and 
spin-spin splitting will decrease, since these all arise from short-range 
effects. Therefore it is to be expected that slopes will become parallel for high 
excitations, where there are many nodes in the radial wave functions. At the 
lowest masses, the tt, 77 and 77' have been omitted from the trajectories because 
of the likehhood that their masses are affected by the instanton interaction. 

For 2+ states, it is rather obvious from Fig. 1 that states lie systematically 
higher than ^P2- A sphtting of ~ 80 MeV is observed. This is natural. For the 
F-states, there is a centrifugal barrier at short range. The need to overcome 
this barrier by attraction from the confining potential pushes the masses of 
F-states up. According to this argument, masses of D states should lie ~ 40 
MeV below F-states. The mass of the best established D-state, p3(1982) is 
indeed 36 ± 15 MeV below that of /4(2018). 

The ccntroids of ^^2, ^-^3 and states with / = 0, weighted by their multi- 
phcities, are at 2292 ± 15 and 2024 ± 8 MeV. The ^F^ states with / = are at 
2275 ±25 MeV and 2025 ±20 MeV. Within the sizable errors, these are consis- 
tent with the absence of spin-spin splitting. The well estabhshed 7^2(2267) hes 
25 ±21 MeV below the centroid of the triplet F-states; and the somewhat less 
well determined 772(2005) lies 19 ± 18 MeV below the centroid of the F-states. 
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However, there is the possibihty that the 772(2005) may be moved upwards by 
the well known level repulsion by 772(1880). In any case, these determinations 
are not as definitive as the difference measured by p3(1982). 

One would then expect G-states to lie ~ 40 MeV above F-states. The situation 

is not clear here because most G-states are located with rather large errors. 
In the absence of tensor and spin-orbit splitting, ^F^, ^F^ and ^F2 states are 
degenerate in mass. 

Tensor and spin-orbit sphtting may be assessed from the relations [79] 



For both multiplcts with I = 0, C = +1, states lie highest in mass, 
requiring significant tensor splitting. This splitting is in the sense predicted 
by one-gluon exchange. Using a linear confining potential plus one-gluon ex- 
change, Godfrey and Isgur [80] predict SMt — 9 MeV for the lower multiplet. 
From Table 5, one finds 5Mt = 20 ± 4.5 MeV for the lower multiplet and 
7.1 ± 7.3 MeV for the upper one; errors take into account observed correla- 
tions in fitted masses. Because the same centrifugal barrier is used for ^^4, 
^Fa and ^F2 states, there is almost no correlation of AMt or AMls with the 
radius of the barrier. The observed spin-orbit sphttings are SM^s = 2.4 ± 4.4 
MeV for the lower multiplet and —6.3 ± 5.1 MeV for the upper one; these 
average approximately to zero. 

Glozman [81,82] has associated the small spin splitting at high masses with 
Chiral Symmetry restoration and the string picture of hadrons. He argues 
that this leads to parity doubling high in the spectrum and cites a number of 
examples. 

It is well known that Regge trajectories arise naturally from the string pic- 
ture of hadrons: quarks tied together by a flux-tube of gluons. The argument, 
originally due to Nambu, is in outline as follows [83]. Suppose the mass origi- 
nates from a uniform fiux tube joining two gluons. Suppose also that for light 
quarks the system rotates and the quarks move at the speed of light at radius 
R. Then the mass associated with the flux tube is 



AMls = [-20M{^F2) - 1M{^F^) + 27M(=^F4)]/54, 
AMt = [-4M(3F2) + 7M{^F^) - ^M{^Fi)]/U. 



(21) 
(22) 




(23) 



The angular momentum is given by 




(24) 
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Fig. 56. Regge trajectories. 
Then J = M^/2T:a, where a is the density per unit length. 

Fig. 56 shows Regge trajectories constucted from mesons quoted in this review. 
On Regge trajectories, is plotted against J. Here the centrifugal barrier 
comes into play, so the slopes of these trajectories will be slightly larger than 
for those plotted against radial excitation. It is found to be 1.167 GeV^. 

It is not clear how to assign 1~ states to these trajectories. The p(1450) does 
not seem to fit naturally as the first radial excitation of the p(770). It has been 
suggested by two groups that the mass of the p' should be lower, ~ 1.29 GeV 
[84,85]; this has been a well known controversy for some years. The lower mass 
would fit much better to a trajectory through the well defined p(1982). Then 
p(1700) and p(2000) fit on to the next two trajectories. However, the p(1700) 
is approximately degenerate with ^D^ and is often taken to be the ^Di state. 
The p(2000) is found experimentally to have a significant ^Di component, but 
could be a mixed state. The 1 states present a dilemma. 

Better data are needed for 1 states. A way forward would be to study 
diffractive dissociation of photons to distinctive channels like 6i(1235)7r. This 
decay can go via L = or 2. It is likely that matrix elements for these two 
decays will be distinctively different from initial states ^Di and ^Si, because 
of the D-wave tied up in the ^Di state. At higher energies, decays to pa 
and p/2(1270) are likely to throw light on the D-wave content of the initial 
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state. If the luxury of a linearly polarised photon beam is available, D-wave 
amplitudes will appear directly via azimuthal dependence of the form cos(f) 
and cos 20 with repect to the plane of polarisation. This arises through the 
partial wave decomposition 

|1, 1 >= ^-12,2 >^ |1,-1 >, Y^|2, 1 >L |1,0 >, +y Y5|2,0 >L |1,1 >. 

A linearly polarised photon is a superposition of initial states |1,1 > and 
|1, —1 > and dependent terms arise in the differential cross section from the 
real part of the interference between these components. 



10 Glueballs 



The latest Lattice QCD calculations of Morningstar and Peardon [86] predict 
a sequence of glueballs beginning with J^^ = 0++ at about 1600 MeV, 2++ 
around 2100 MeV, then 0~+ and a second 0+"*" close by in mass. There are 
'extra' states which do not fit naturally into qq spectroscopy at 1500, 1980, 
2190 and 2100 MeV with just these quantum numbers. All except the /2(1980) 
appear prominently in J/^ radiative decays, the classic place to look for glue- 
balls; and there is tentative evidence for the /2(1980) too. 

Lattice QCD calculations suffer from a scale error of ~ (10 — 15)%, so one 
should be prepared for some overall shift in mass scale. Also the Lattice calcu- 
lations are presently done in the quenched approximation, where coupling to 
qq states is omitted. However, mass ratios are predicted with smaller errors. 
Table 10 compares mass ratios of the four candidate glueballs with Lattice 
predictions. There is remarkable agreement, suggesting this is not an acciden- 
tal coincidence. This coincidence was pointed out by Bugg, Peardon and Zou 
[87] . Here the properties of the candidate glueballs will be reviewed in detail. 



Ratio 


Prediction 


Experiment 


M(2++)/M(0++) 


1.39(4) 


L32(3) 


M(0-+)/M(0++) 


1.50(4) 


L46(3) 


M(0*++)/M(0++) 


L54(ll) 


L40(2) 


M(0-+)/M(2++) 


L081(12) 


1.043(36) 



Table 10 

Glueball mass ratios predicted by Morningstar and Peardon [86] compared with 
experimental candidates described in the text; errors are in parentheses. 

What is clear immediately is that the candidate states do not decay flavour 
blind, except perhaps for the 0"+ state. There is obviously mixing with nearby 
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qq states, as one would expect from the fact that gluon and quark jets develop 
readily at high energies. At present, this mixing is still being unravelled. 

10.1 /o(1500) and the qq nonet 

The evidence for the three 0"^ states /o(1370), /o(1500) and ao(1450) appeared 
in four papers of Anisovich et al. early in 1994 [88-91]. In the first of these, 
it was already realised that /o(1370) and /o(1500) lie too close together for 
them both to be non-strange qq states. The /o(1370) was associated with a 
nonet containing also A'o(1430) and ao(1450); the glucball identification of 
/o(1500) was pointed out. Ref. [90] gave a detailed exposition of the N/D 
technique used in the analysis and Ref. [91] gave extensive details of results. 
These results were confirmed 18 months later with higher statistics and using 
K- matrix techniques by Amsler et al. [92]. 

With the benefit of hindsight, it is clear that /o(1370) was observed as the 
S or e(1300) meson in earlier experiments on tttt — > KK [93-96], where a 
small but definite peak occurs in the S-wave at 1300 MeV. It also appeared 
in Att in the work of Gaspero [97] on pn 27r"'"37r~. Establishing masses 
and absolute branching ratios of nonet members is important in isolating the 
glueball component. The /o(1370) decays weakly to nn, rirj and KK; decays to 
Att have a branching ratio at least a factor 5 larger than to tttt, as demonstrated 
in the analysis of Svr data from Obelix [98] and data on tjAtt [99] and Svr [100- 
102] from Crystal Barrel. In view of the dominant An decay, it is unlikely to 
be an ss state. 

An important point in analysing data is that the peak observed in An is shifted 
up by ~ 75 MeV with respect to that in 27r, because of the rapidly increase 




M (GeV) M (GeV) 



Fig. 57. (a) The shape of /o(1370): in the 27r channel (chain curve), in the Att channel 
(dashed); the full curve shows a Breit-Wigner resonance of constant width, (b) m(s). 
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of 47r phase space with mass. This is illustrated in Fig. 57(a) [99]. The chain 
curve shows the expected 27i cross section and the dashed curve shows that 
expected for An; all curves are normalised to 1 at their peaks. Also there has 
been no full analysis yet including the term m{s) of the full Breit-Wigner 
amphtude of eqn. (4) of Section 4.1; it arises through dispersive effects from 
the s-dependence of r(47r) and is shown on Fig. 57(b). It has a large effect 
and consequently all claims about the mass and width should be treated with 
some reserve. The /o(1370) does not appear significantly in the Cern-Munich 
data for tttt elastic scattering [103,104], where its intensity depends on r|^. It 
does appears clearly in the new BES II data on J/\t' ^ ^tttt. Fig. 26, Section 
5.5; the mass found there is 1335 ± 27 MeV and the width 265 ± 60 MeV [37]. 

With the benefit of hindsight, it is also clear that the /o(1500) was observed 
by Kalogeropoulos and collaborators [105-7]. It was also observed as ^(1590) 
by GAMS [108-112], but shifted in mass through interference effects; they 
appreciated its significance and claimed it as a ghicball because of its rji] 
decay, regarded then as a touchstone for gluonic content. 

Since its discovery, the /o(1500) has been identified as a An peak in a re-analysis 
of Mark III data on J/^ ^An [39]. The /o(1770) and /o(2100) appear as 
two further peaks. The data have been shown earlier in Fig. 25, Section 5.5, 
together with E760 data which show a striking resemblance. In the original 
Mark III and DM2 analyses [113,114], the peaks were not separated from a 
large, slowly varying pp component with — 0~. The analysis of J/^ — > ^An 
has been confirmed by BES I [115]. Fig. 58, taken from this analysis, shows 
the slowly varying 0~ signal in (a), and 0"*" and components in (b) and (c); 
these are determined by the spin-parity analysis and are not conspicuous in 
the raw data shown as points with errors. The latest BES II data [116] have 
7 times higher statistics and confirm the analysis, improving on all individual 
components. 

From the early analysis of J/^f radiative decays, there was some confusion 
whether the fj at 1710 MeV in KK has J = or 2. This has now been 
estabhshed as J = by a re-analysis of those data by Dunwoodie [117]. It has 
also been confirmed by the Omega group [118] and from new BES II data on 
J/^ ^-fKK [119]. 

The ao(1450) has recently been confirmed in decays to cop in Crystal Barrel 
data on pp — > a;7r"'"7r~7r° at rest [120]. Decays to cop are a factor 10.7 ± 2.7 
times stronger than to rjn, where is was originally discovered. The weak decay 
to r^TT explains why it has not been observed there in E852 and VES data. Fig. 
59(a) shows the percentage contribution of ao(1450) as its mass is scanned; 
(b) shows the line shape in different channels; the peak appears quite narrow 
in the up channel because of the unavoidable effect of the rapidly rising r(pci;) 
V. mass. The percentage contribution fitted to data is shown in (c) v. width 
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Fig. 58. Components from the BES I analysis of J/^ — > 7(47r). 

when a Breit-Wigner amphtude of constant width is used as a test; (d) shows 
the Argand diagram in the r/vr channel. 

In conclusion, there is available a standard nonet made up of /o(1370), ao(1450), 
ii'o(1430) and /o(1710), similar to known and ^Pi nonets. There is some 
doubt within ~ 70 MeV of the precise mass of /o(1370): two-body data tend 
to favour masses of 1300-1335 MeV and 4-body data favour higher masses 
close to 1400 MeV. The mass of ao(1450) may be affected by its strong cou- 
pling to the up threshold. So it would be a mistake to make too much of mass 
differences amongst the lowest three members of this nonet. The /o(1710) is 
rather further in mass from /o(1370) than /2(1525) is from /2(1270), but the 
/o(1370) and /o(1710) could be pushed apart by the usual level repulsion due 
to an extra state /o(1500) between them. 



10.2 Mixing of the glueball with qq states 



There have been many studies of the mixing, based on both the observed 
masses and on branching ratios of states between vrvr, rji], rji]' and KK. The 
much stronger coupling of /o(1500) to vrvr than to KK implies that it is not 
primarily an ss state. Amsler and Close [121] consider the mixing of nn and 
ss components into the glueball. If the composition of the mixed state \G > 
is written 

,^ ^ _ \Go > +^[V2\nn > +uj\ss >] 
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Fig. 59. (a) Percentage contribution of ao(1450) v. fitted mass M using a Flatte 
formula including r(r/7r), T{puj) and m(s). (b) The line-shapes for the Flatte am- 
plitude itself (full curve), its contribution to cop (dashed) and its contribution to 
pp — > {ujp)tt (dotted), (c) Percentage contribution of ao(1450) v. fitted width for 
a Breit-Wigner amplitude of constant width, (d) the Argand diagram for the rjir 
channel; numerical values of mass are shown in GeV. 

then lowest order perturbation theory gives 

^^<dd\V\Go> /{Ea^-Enn) (26) 
u = {Eg„ - Er,n/{EG„ - E,,), (27) 

where E are masses of unmixed states, Gg refers to the unmixed glueball and 
V is the interaction responsible for mixing. If Eq^ lies midway between the 
masses of unmixed nn and ss states, \G > will acquire extra dd component, 
losing ss component. The state G then moves towards a flavour octet, in 
agreement with observation for /o(1500). If the bare glueball were to lie close 
to /o(1710), it would acquire a large ss component and move towards a flavour 
singlet. One does not see two neighbouring ss states. This suggests that the 
mass of Go > lies close to 1500 MeV. 

Close and Kirk [122] use observed masses and WA102 branching ratios and 
arrive at wave functions as follows: 

|/o(1710) > = 0.39\G > +0.91|ss > +0.14|nn > (28) 
|/o(1500) > = -Om\G > +0.37|ss > -0.62|nn > (29) 
|/o(1370) > = O.eOlG > -0.13|ss > -0.79|nn > . (30) 
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Strohmeier-Presicek et al. [123] arrive at somewhat different numerical coeffi- 
cients, but tlie overall picture is similar. All authors neglect mixing between 
/o(1710) and the nearby /o(1770). Amsler points out that the observed upper 
limit for /o(1500) — >• 77 is inconsistent with a pure qq state [124]. 

In a series of papers, Anisovich and Sarantsev examine a somewhat different 
scheme [125], using the K- matrix. They take the view that K-matrix poles 
represent 'bare', unmixed states. The coupling constants of these poles to tttt 
etc. then 'dress' these bare states and allow for the effects of meson decays. 
The virtue of their approach is that they fit a very large volume of data on tttt 
elastic scattering, tttt — >■ KK and 7777, plus many channels of pp annihilation 
at rest. There is also the virtue that bare nn, ss and \Go > states are fitted 
directly to data. They arrive at five K-matrix poles from 650 to 1830 MeV and 
identify the lowest nonet with T-matrix poles as /o(980), ao(980), ^^0(1430) 
and /o(1370), all strongly mixed with a very broad glueball with mass 1200- 
1600 MeV. The /o(1500), ao(1450) and /o(1710) then fall into the first radial 
excitation. 

There can be no doubt that the amplitudes fitted to data reproduce the same 
narrow T-channel resonances and the same branching ratios as are found in 
other analyses. The difference lies in how to interpret the broad component 
which appears so conspicuously in the tttt S-wave from threshold to 2 GeV and 
can be built into other channels also as a broad underlying background [16]. 
I have had extensive and friendly discussions with Anisovich and Sarantsev 
about this scheme. It has the awkward feature that the lowest K-matrix pole is 
largely ss; an obvious question is whether this could be replaced by the a pole 
discussed in the next section. Another point is that the broad structure at high 
mass may arise from the opening of the strong Att threshold and associated 
dispersive effects. This requires fitting An cross sections to data. It is rather 
important to establish whether their scheme is unique. An alternative is to fit 
(i) a, K, /o(980) and ao(980) as the lowest nonet (not necessarily qq), (ii) the 
conventional nonet of /o(1370), ao(1450), Ko{lA30) and /o(1770) mixed with 
/o(1500), plus (iii) /o(1770) as a member of the next nonet. The question is 
whether this rearrangement of the broad background fails to fit the data, and 
if so how. 

In all schemes, there is agreement that the data favour the presence of a 
glueball in the mass range close to 1500 MeV. Given this and the Lattice 
QCD calculations, the question is: where are the next three glueballs? This 
question is addressed in the following. 
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Fig. 60. The mechanism of central production. 
10.3 Features of Central Production Data 

In central production, two protons scatter through a small amgle with small 
momentum transfer. The standard hypothesis for the mechanism of the re- 
action is shown in Fig. 60. Two gluons, or a ladder of gluons making the 
Pomeron, undergo a collision close to the centre of mass of the system and 
produce a meson R. There are extensive data from the Omega and WA102 
collaborations at CERN for many final states. 

Close and Kirk found a mechanism for filtering out well-known qq states by 
selecting coUisions where the two exchanged gluons have a small difference 
in transverse momentum, APt [126]. Well known qq states such as /i(1285), 
/2(1270) and /o(1370) are P-states and have a zero at the origin of their wave 
functions. The idea is that they will be suppressed close to zero transverse 
momentum in the collision forming resonance R. It is observed experimentally 
that this does indeed suppress well known qq states, leaving exotic states such 
as /o(1500). 

The idea developed further by examining the angular dependence of the col- 
hsion. Treating it in analogy to e'^e~ — > R, Close, Kirk and Schuler studied 
the azimuthal dependence of the production process [127-129]. They define (f) 
as the azimuthal angle between the transverse momenta of the two exchanged 
gluons (reconstructed from the kinematics of the outgoing protons). From the 
model, the azimuthal dependence is predicted uniquely for production of 0^, 
1+ and 2~ resonances R. The data agree closely with predictions, vindicat- 
ing the analogy between e~^e~ interactions where photons are exchanged and 
central production where gluons are exchanged. 

If R has quantum numbers 0+ and 2+, the azimuthal dependence is not unique. 
It depends on two cross-sections which they denote as TT and LL, standing for 
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Fig. 61. The (j) dependence for central production of /o(1370), /o(1500), /2(1270) 
and /2(1950). 
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Fig. 62. WA102 data for central production of tt+tt tt+tt : (a) total mass spectrum, 
(b) 1++, (c) 0++ pp, 0+ (Tcr, (e) 2-+ and /2(1270)7r7r compared with the fits. 

transverse or longitudinally polarised gluons. For 0+ and 2+ mesons, results are 
shown in Fig. 61 and are truly remarkable. The well known qq states /o(1370) 
and /2(1270) have one distinctive form of angular dependence; /o(1500) and a 
broad 2+ state at ~ 1950 MeV have a quite contrary azimuthal dependence. 
Although the origin of this difference is not understood yet, it suggests some- 
thing different about /o(1500) and /2(1980). If /o(1500) has a large glueball 
content, it suggests that perhaps /2(1950) does also. 
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The broad /2(1950) in central production data of 47r was first reported by the 
Omega collaboration in 1995 [130] with a width of 390 ± 60 MeV. It may be 
the same object as reported earlier in K*K* by LASS [131] with a width of 
250 ± 50 MeV. The WA102 collaboration produced higher statistics in 1997 
[132] and again in 2000 [133] in 47r°, 7r+7r-7r°7r° and 7r+7r-7r+7r-. The latest data 
are shown in Fig. 62. Decays are consistent with being purely to /2(1270)(T, 
where the a has parameters close to those of the a pole discussed in the next 
Section. The An peak has M = 1995 ± 22 MeV, T = 436 ± 42 MeV for data 
containing 4 charged pions and M = 1980 ± 22 MeV, T = 520 ± 50 MeV for 2 
charged pions. Branching ratios between the three sets of data are accurarely 
consistent with 7 = 0. Note in passing that Fig. 62(e) confirms the 7^2(1645) 
and 772(1860) discussed earlier in Section 5.3. Also note that the O"*""*" signal in 
(d) is at 1500 MeV and narrow. In (c) there is a broader 0"*"+ signal peaking 
at 1450 MeV and coming from interference between /o(1370) and /o(1500); 
the upper peak in this figure is evidence for the 0+ state called /o(2020) in 
Section 5. 

Meanwhile Crystal Barrel observed a very similar broad 2+ peak in tjt] with 
M = 1980 ± 50 MeV, F = 500 ± 100 MeV, values very close to WA102. The 
observation was made in pp — > rjrjn^ at momenta 1350 to 1940 MeV/c. The rfri 
decay angular distribution is shown in intervals of rji] mass in Fig. 63. In the 
lowest interval, 1.45-1.55 GeV, the process is dominated by a strong /o(1500) 
signal and the decay angular distribution is almost flat. At higher masses one 
sees the gradual onset of a higher partial wave which is fitted as purely 2+; no 
4+ component is observed. The gradual onset of the 2+ amphtude requires a 
very broad structure inconsistent with /2(1934). 

A curious feature of the data was reported in this publication, but not un- 
derstood at the time. The broad /2(1950) is produced with orbital angular 
momentum L = 1 in the final state and with component of spin ±1 along 
the beam direction. Later, when the combined analysis of/ = l,C = +l 
was carried out, the explanation of this result was uncovered. The /2(1950) is 
produced purely from the decay of the a2(2255) in the rjrjn'^ final state. The 
value of the ratio r2 F2/^P2 for that resonance fits as r2 = — 1.9 ± 0.4. 
Clebsch-Gordan coefficients for coupling of pp to ^^2 and are such that 
the initial state will have only helicity if r2 = — y^7/2. This helicity state de- 
cays purely to final states with helicity ±1, again because of Clebsch-Gordan 
coefficients. This explains the curious observation of the helicity of /2(1950). 
However, it also firmly estabhshes that 02(2255) is dominantly ^^2. 

The broad /2(1950) makes a large improvement to in the combined analysis 
of J = 0, C = +1 data [3] in the rjrj channel. Without it, there are awkward 
interferences between the two and ^^2 interferences, so as to reproduce 
not only the narrow states but generate also a broad background. Adding the 
broad background explicitly has the effect of decouphng the narrow states to a 
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large extent. The broad component fits with M = 2010±25 MeV, F = 495±35 
MeV. It appears also in tttt with gf^^/g^o^o = 0.72±0.06. This is to be compared 
with the prediction 0.41 for nn and 1 for a glueball. It leads to a fiavour mixing 
angle $ = (23.6 dz 3.5)°, compared with the value 35.6° for a glueball. Some 
mixing between a glueball and neighbouring qq states is plausible. 

There does however remain one anomaly which needs further work. The WA102 
collaboration has studied central production of rjfj. There they do not find the 
broad /2(1980) or /o(2100) which are so conspicuous in Crystal Barrel data 
[134]. Instead they observe an /2(2150) peak with M = 2151 ± 16 MeV and 
width r = 280 ± 70 MeV; but on close examination the data actually peak 
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Fig. 64. Etkin et al. data on 0^: (a) acceptance corrected mass spectrum, (b) inten- 
sities of partial waves, (c) relative phases; curves show the fit to three resonances. 

at 2060 MeV because of interference with the S-wave. The actual moments 
tiM fitted to the data are quite noisy, see their Fig. 3. This discrepancy is not 
understood at present. 

There is one small comment on the mass. The PDF quotes an average mass of 
1934 ± 12 MeV. This average is strongly influenced by the low value reported 
by Amsler et al. [49]: M = 1867 ±46 MeV. The data used in that publication 
were at 900 MeV/c where the upper limit available for 7777 mass is 1914 MeV 
and the broad 2"*" signal is much weaker than in the momentum range 1350 
to 1940 MeV/c. If one instead averages the results of Anisovich et al. and the 
post-1995 Central Production data, the average mass is 1972 ± 12 MeV. 

10.4 Tin — > (j)4> 

Etkin et al. [135] observe J^*^ = 2++ peaks in tttt (f)(1) at 2150 and 2340 
MeV, Fig. 64. This process violates the OZI rule. The overall observed cross 
section is a factor 100 stronger than can be explained readily by OZI conserving 
processes. An obvious suggestion is that the peaks are due to SU (3) singlet 
states which mix strongly with a broad 2+ glueball. The narrow states can be 
produced via an nn component and decay via the ss component. An /2(2150) 
is listed by the PDG, corresponding to observation by WA102 in 7777 [136] and 
KK [117] and by GAMS in tjt] [137]. However, remaining entries listed by the 
PDG were quoted as the /o(2100) in the original pubhcations. 

An ss state is expected ~ 250 MeV above /2(1934) and could well explain 
the /2(2150). It would fit in with the S-wave decay to 00 of the 2150 MeV 
peak. Then a ss state is expected ~ 100 MeV higher. Etkin et al. claim 
two D-wave (jxj) states at 2297 ± 28 MeV (with P = 149 ± 41 MeV) and at 
2339 ± 55 MeV (P = 319tg9) MeV. An obvious question is whether there is 
really one D-wave state and an interfering broad 2+ background from /2(1950). 
A re-analysis of the data is required to answer this question. 

Palano [138] has presented JETSET data taken at LEAR on pp (jxf), Fig. 
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Fig. 65. JETSET data for pp cfxp. 

65. A peak is observed in the DO wave (J = 2, 5* = 0, L = 2). It is fitted with 
M = 2231 ±2 MeV, F = 70± 10 MeV; however there is noise in the signal and 
Palano remarks that results could be compatible with the Etkin et al. peak 
at 2297 MeV. This is because of interference with a broader stucture in the 
D2 wave, J = 2,5' = 2,L = 2, in the first panel. Phase motion is observed 
between the DO and D2 components. 

10.5 Evidence for the 0~ glueball 

This is an entertaining story of unitarity and analyticity at work. It has ap- 
peared in two publications by Bugg, Dong and Zou [139,140]. However, even in 
these publications the full ramifications of the story were not fully appreciated 
and a further twist needs to be given here. 

Radiative J/^ decays are a classic hunting ground for glueballs. The J/^ has 
process which conserves C-parity is the radiation of one photon 
and two ghions. The two gluons will couple preferentially to glueballs because 
«5(M^) ~ 0.5 at these masses, so the intensity of coupling to qq states is 
suppressed by a large factor. 

All of the radiative decays to rjnn, KKti, in and i^*(890).ft'(890) contain a 
dominant 0^ final state. This is illustrated in Fig. 66 with data from a variety 
of experiments. In (a), from Mark III [141], attention was focussed on the two 
narrow peaks at low mass. What is not displayed is a broad 0^ signal under 
the peaks, covering the whole mass range; this broad component is shown on 
Fig.58(a). Other panels of Fig. 66 show BES I data [115,142,143]. In aU three 
cases there is a broad 0~ peak. Wermes made the first attempt to fit this broad 
signal [144]. 

Bugg, Dong and Zou [140] showed that it may be fitted by a single resonance. 
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Fig. 66. (a) Mass spectrum for J/^ radiative decays to (a) (rjir'^Tr ) after selecting 
ao(980) [141], (b) K*K* [142] and (c) K^K^tt^ [143]. 
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Fig. 67. The line-shape of the broad component after dividing out the E'^ factor. 

as will be explained below. However, their paper failed to point out one im- 
portant feature. In the production process J/\l/ — 70^, there is one unit of 
orbital angular momentum for production. This leads to an factor in the 
intensity, where E.y is the photon energy. This factor inflates the lower side of 
the resonance. Fig. 67 shows the line-shape after dividing out the factor E^. It 
peaks at 2260 MeV with a conventional full width at half maximum of ~ 350 
MeV. The shoulder on the lower side will be explained below. 

A single Breit-Wigner form is used to fit simultaneously all decay channels i: 



fr. 



(31) 



M^-s- m{s) - iM ^i{s) ' 
where Aj are coupling constant for each channel; rj(s) are taken to be pro- 
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Fig. 68. Widths Ti{s) for A) j^vrvr, B) pp, C) uu, D) K*K*, E) 00, F) kK and G) 
total; C is shown dashed for clarity. 

portional to the phase space available for decay to each channel, multiplied 
by a form factor exp(— i?^p^/6) with R = 0.8 fm. These widths are shown in 
Fig. 67 after normalisation to cross sections in Fig. 66. The function m(s) is 
an s-dependent dispersive correction to the mass: 



The total width Ttot tends to zero at high mass sufficiently rapidly that the 
dispersion integral converges well with a subtraction at s = M^. 

On Figs. 69, the magnitude of each Fj is fitted freely. The fit gives coupling 
constants squared for pp, ujuj, K*K* and (fxf) remarkably close to the ratio 
3:1:4:1 expected for fiavour-blind decays. However, the available phase space 
has quite a significant impact: obviously, the 00 channel is forbidden below its 
threshold, and kaonic channels are significantly suppressed. The fitted mass is 
2190 ± 50 MeV. 

Fig. 69 shows the actual fit to data, where the E'^ factor is included. This figure 
shows the essential features of the fit. At low masses, the only open channel is 
T^TTTT. Then strong channels pp and ujuj open at ~ 1550 MeV. Their widths in 
the denominator of the Breit-Wigner amplitude cut off the ?77r7r channel rather 
sharply, as one see in Fig. 69 (a). This cut-off is responsible for the shoulder 
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{s) = is-M') I 



(32) 
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Fig. 69. The fit to the broad component of the signal in (a) r/vrvr, (b) pp, (c) ujlo^ 
(d) K*K*, (f) (/)(/) and (f) the sum of all channels [140]. 

in the total cross section of Fig. 67. Likewise, branching ratios to pp and uouj 
are cut off at 1750 MeV by the opening of the strong K*K* channeL The 
resulting line-shape is rather difficult to interpret. Perhaps one should take a 
somewhat larger full width than 350 MeV in order to allow for the shoulder at 
low mass. The fitted mass 2190 ± 50 MeV is somewhat below the peak value 
because of the strong effect of the K*K* width. 



10.6 Sum Rules for J/\E' radiative decays 



Branching fractions for production of glueballs in J/\l' l{gg) will now be 
considered. Close, Farrar and Li [145] make predictions for branching fractions 
for individual of gg; results depend on masses and widths of resonances. 
Using experimental values for these. Table 11 [87] makes a comparison of their 
predictions with experiment. 
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State 


Prediction 


Observation 


/o(1500) 


1.2 X 10" 


-3 


(1.03 ±0.14) X 10-^ 


r?(2190) 


1.4 X 10" 


-2 


(1.9 ±0.3) X 10-2 


/o(2105) 


3.4 X 10" 


-3 


(6.8 ± 1.8) X lO"'' 


/2(1980) 


2.3 X 10" 


-2 


(2.6 ±0.6) X 10-3 



Table 11 

Branching fractions of glueball candidates in J/^ radiative decays, compared with 
predictions of Close, Farrar and Li [144]. 

The branching ratio for /o(1500) is close to prediction for a glueball and far 
above that predicted for qq. The production of 77(2190) in pp, uu, (f)(f), K*K*, 
Tjinr and KKn [14,28] is slightly above the prediction. The /o(2105) has so far 
been observed in J/\I' radiative decays only to aa [14]. For this mode alone, 
the observed production is only 20% of the prediction; decays to r]r] and tttt 
will also contribute, but their decay branching ratios relative to Att are not 
presently known. 

The prediction for the branching fraction of the 2+ glueball is large if the 
width is taken to be the 500 MeV fitted to /2(1980). Observed decays to aa 
and_/2(1270)(7 account for (10 ±0.7 ±3.6) x lO"'' of J/* radiative decays and 
K*K* decays a further (7 ± 1 ± 2) x 10"^. If one assumes flavour- bhndness 
for vector-vector final states, the vector-vector contribution increases to (16 ± 
2 ± 4.5) X 10""^. The total of 2.6 x 10"^ is stiU a factor 9 less than predicted 
for a glueball. This is presently a problem in identifying /2(1980) with the 2"*" 
glueball. The total prediction for glueball production in J/^ radiative decays 
is 4.2%, compared with ~ 6% for all radiative decays. As yet, only half of the 
products of J/^ radiative decays have been assigned to specific J^. Data on 
radiative decays to rjrj and rjfjTn: would be particularly valuable. 

Longacre [146] has made an interesting suggestion about the missing 2+ com- 
ponent in J/^ radiative decays. He suggests that the glueball may be a ring, 
like a doughnut, with a hole in the middle. If so, it has a D-state wave fimc- 
tion. There would then be a centrifugal barrier depressing its production rate 
in J/^ radiative decays. A better place to look would be in ^" decays, where 
the initial state has a D-state wave function providing a better match to the 
glueball. The channel to study is ^" — > 77(47r), but there are no data at present. 

However, it is not obvious that this conjecture is consistent with central pro- 
duction data. Close and Kirk argue that qq states are suppressed at small 
APt because they have P-state wave functions. A D-state glueball would be 
suppressed even further. Possibly there is some other explanation of the en- 
hancement of /o(1500) and /2(1950) at small APt in central production. 
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11 The light scalcirs 



There has been much argument whether a and k exist; also whether /o(980) 
and ao(980) are part of the ground-state qq family or whether they are 4-quark 
states or KK 'molecules'. Over the last 2 years valuable new information have 
appeared. The experimental situation will be reviewed first. This is followed 
by a brief summary of theoretical approaches. A rough concensus is that the 
light scalars make a nonet separate from qq states, but probably with some 
qq mixed in and probably with significant breaking of SU(3) due to meson 
masses. 

The E791 collaboration has observed a low mass tt+tt^ peak in tt+tt^tt"'" 
[147]; they fit it with a simple Breit-Wigner resonance with M = 478li MeV, 
r = 324l|g ±21 MeV. A low mass peak was observed earlier in DM2 data for 
J/^ iOT^^T^~ [148]. New BES 11 data [149] also contain a similar peak shown 
in Fig. 70(c). The o appears as a conspicuous band at the upper right-hand 
edge of the Dalitz plot in (b). The fitted 0"'' and 2+ components are shown 
in (e) and (f); the cross-hatched area in (c) shows the coherent sum of a and 
and a small /o(980) contribution. 

Why does the a not appear the same way in tttt elastic scattering? The an- 
swer is discussed in a recent paper which refers to both a and k [150]. In tttt 
elastic scattering, there is a well known Adler- Weinberg zero just below the 
elastic threshold [151] at sa — ''^^/2. This zero is today a central feature of 
Chiral Perturbation Theory. To a first approximation, it cancels the pole in 
the denominator near threshold, with the result that 6{s) is rather featureless 
at low mass and shows little direct evidence for the pole. 

This Adler zero is absent in production processes. The production amplitude 
may be taken to be 

where F is s-dependent. For tttt elastic scattering, the amplitude must lie on 
the unitarity circle and the corresponding amplitude is 

r(7r7r ^ tttt) = ^^^^Lr.. ^ ■ (34) 

This contains an extra factor MY{s) in the numerator. The denominator is 
common to all tttt reactions. If a Breit-Wigner amplitude describes the tttt 
elastic amplitude, numerator and denominator must be self-consistent. 

For elastic scattering, there are other small contributions from higher reso- 
nances /o(980), /o(1370), /o(1500), etc. They are accomodated by adding their 
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Fig. 70. (a) The uj peak; the lower curve shows the quadratic fit to the background, 
(b) the Dahtz plot for ujit^tt^. (c) and (d) projections on to tttt and lott mass; his- 
tograms show the fit and the shaded region the background estimated from sidebins; 
the dashed curve in (d) shows the fitted 61 (1235) signal (two charge combinations), 
(e) and (f) show mass projections of /o and /2 contributions to 7r"'"7r~ from the fit; 
in (e), the shaded area shows the a contribution alone, and the full histogram the 
coherent sum of a and /o(980). 
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Fig. 71. (a) Phase shifts for tttt elastic scattering; dark points are phase shifts of 
Hyams et al. [104] above 450 MeV and from Ke4 [152] data below 400 MeV; open 
circles show data [153,154]; the curve shows the fit. (b)The trajectory in the 
complex s-plane where the real part of the amplitude goes through zero. 

phases to that of the a; this prescription satisfies elastic unitarity. In principle 
there could be a further background amplitude. The question is whether the 
simple prescription of eqns. (34) and (35) fits the data. In practice it does. 

The rate for a transition is proportional to \F\^, where F is the matrix element 
for the transition. The width of a resonance is proportional to this transition 
rate. It is therefore logical to include the Adler zero explicitly into the width. 
The idea behind this is that the matrix element for tttt elastic scattering goes to 
zero for massless pions as /c — > 0. In production reactions such as J/^' — > cua, 
the pions are no longer 'soft' and no Adler zero is required. 

The parametrisation used in Ref. [150] is 

\ f \ ^ ^ ff \ -(s-M^) 

r(s) = p(g) ^2 — J — ■ (35) 

Here sa — "1-^/2, the position of the Adleer zero. Empirically, the best form for 
f{s) is bi + 62S. A combined fit is made to BES II data, CERN- Munich phase 
shifts for TT^TT^ — TT"?!"*" [104], the most recent Ke4 data [152], CP violation 
in decays [11] and data for n~7i~^ — > n^n^ [153,154]. The fit to elastic data 
is shown by the curve in Fig. 71(a). There are some systematic discrepancies 
above 740 MeV between phase shifts of Hyams et al. (filled circles) and n^n^ 
data (open circles) and the fit runs midway between them. Fitted parameters 
are M = 0.9254, A = 1.082, bi = 0.5843, 62 = 1.6663, all in units of GeV. 

A consistent fit is found to both elastic and production data. As a check, the 
BES data alone have also been fitted with (i) F = Fop(s)/p(M^), where p(s) 
is 2-body phase space, (ii) a Breit-Wigner amplitude with constant width. 
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All three pole positions shown in Table 12 are in close agreement. The first 

entry is the most reliable since it fits all data. The square of the modulus of the 
amplitude in cqn. (34) peaks at 460 MeV, slightly below the peak in Fig. 70(c) 
because of the phase space for ujtttt at low tth mass. The pole is in reasonable 
agreement with the result of Colangelo et al.: M = (490 ± 30) - i(290 ± 20) 
MeV [155]. Their determination has the virtue of fitting the left-hand cut 
consistently with the right-hand cut, using the Roy equations. However, they 
do not fit data above 800 MeV, so the effect of /o(980) is omitted. Markushin 
and Locher [156] show a table of earlier determinations which cluster around 
the values in Table 12, with quite a large scatter. 



Fit 


Form 


Pole position (MeV) 


A 


eqns. (8)-(10) 


533 ± 25 - z(249 ± 25) 


B 


F =cx p{s) 


569 ± 40 - (274 ± 25) 


C 


F = constant 


537 ± 40 - i{268 ± 50) 



Table 12 

Pole positions of the a for the three fits. 

A popular misconception is that the pole needs to be close to the mass M 
where the phase shift passes through 90° on the real axis. That is true only 
when r is constant. Because F varies with complex s in an explicit way, the 
phase of the amplitude varies as one goes off the real s axis to imaginary 
values. Fig. 71(b) shows the position on the complex plane where the phase 
goes through 90°; there is a rapid variation with Im m, related directly to the 
Adler zero in the amplitude near s = 0. 

In BES data, the phase of the a amplitude is determined by strong interfer- 
ences with 6i(1235). There arc two 6i(1235) bands on the Dalitz plot of Fig. 
70(b) and together they account for 41% of the intensity. The a accounts for 
19% and /2(1270) for most of the rest. An independent check is made on the 
phase of the a amplitude by dividing up the tttt mass range into bins 100 MeV 
wide from 300 to 1000 MeV. In each bin, the magnitude and phase of the a 
amplitude are fitted freely, re-optimising all other components. 

Results for the phase are shown on Fig. 72. Curves show phases from the three 

parametrisations described above. [A detail is that zero phase is taken from fit 
A at the tttt threshold. Strong interactions introduce an overall phase for all 
formulae and the dashed and dotted curves have been adjusted vertically to 
agree with the full curve at 550 MeV; only the mass dependence of the phase is 
meaningful.] There is only small discrimination between the three forms. The 
agreement of phases with the curves demonstrates the correlation of magnitude 
and phase expected from analyticity. Magnitudes in individual bins fluctuate 
from 0.90 to 1.17; removing those fluctuations reduces errors on phases by a 
factor ~ 1.4, but does not affect the observed agreement significantly. 
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Fig. 72. The phase of the tttt S-wave amphtude. The full curve shows fit A, the 
dashed curve fit B (Brcit-Wigncr amplitude of constant width) and the dotted 
curve fit C (r oc p{s)); points with errors show results fitted to slices of tttt mass 
100 McV wide. 

There is one technicality concerning the form r(s) = rop{s). It leads to a 
false virtual state pole in the mass range M — 200 — 230 MeV, below the 
TTTT threshold [150, 157]. It is therefore not a good choice of fitting function, 
though Table 12 shows it does not have a drastic effect in fitting data. 



11.1 The K 



The LASS data for Kti elastic scattering [158] are fitted with 

m^pis) roexp-(aM); (36) 

M"^ — sa 

sa = Mj^ + m^/2. As for tttt, it is assumed that phases of the k, and iro(1430) 
add. For Xo(1430), a Flatte form is used: 

r(s) = ^2 _^fj 9lPKn{s) + g2PKr,'{s)]. (37) 

There is some freedom in fitting the 7^0(1430), because it rides on top of the k 
'background'. As other authors have found, there is no improvement if decays 
to Kt) are included. 

Fitted parameters for the n axe M = 3.3 GeV, 61 = 24.53 GeV, a = 0.4 GeV^. 
These parameters are highly correlated and one can obtain almost equally 
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good fits over the range M = 2.1 to 4 GeV by re-optimising bi and A. The 

Ko{U30) is fitted with M = 1.535 GeV, gi = 0.361 GeV, gi/gi = 1.0. Its 
pole position is (1433 ± 30 ± 10) - i(181 ± 10 ± 12) MeV; that is close to the 
values fitted by Aston et al: M = (1412 ± 6 MeV, T = 294 ± 23) MeV. [The 
full width r is twice the imaginary part of the pole position.] 
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Fig. 73. (a) Ktt phase shifts. Points are from LASS [158]. The full curve shows the fit 
to LASS data and the dashed curve the k contribution, (b) The fit to the intensity 
of Ktt clastic scattering from LASS data, (c) The intensity of the k amplitude in 
production processes, normalised to 1 at its peak; the full curve is from the fit to 
LASS data and the dashed curve from the fit to E791 data, (d) the k phase fitted 
to BES II data compared with the curve fitted to LASS data. 

The fitted Ktt scattering length is 0.17m^^, very close to Weinberg's prediction 
from current algebra of 0.18m~^. Fits to Kn phase shifts are shown in Fig. 
73(a) and to the intensity of the clastic amplitude in Fig. 73(b). The pole 
position of the k is at M = (722 ± 60) - «(386 ± 50) MeV. The fitted value of 
M for the k is very high; it is essential that it lies above the i^o(1430) since the 
phase shift does not reach 90° before that resonance. As shown on Fig. 73(d), 
the fitted k, phase reaches ~ 70° at 1600 MeV. It is possible that the phase 
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for real s does not actually reach 90°. A similar case arises for the well-known 
NN ^So virtual state, where the phase shift rises to 60° and then falls again. 

The full curve in Fig. 73(c) shows the intensity of the k, amphtude as it will 
appear in production processes; the curve is normalised to 1 at the peak at 
700 MeV. In production reactions such as — > ir'^K, the observed intensity 
will be multiplied by the available three-body phase space. 

If the Adler zero is removed from the parametrisation of the Ktt S-wave and 
from the ii'o(1430), the k, pole disappears. So the existence of the k, is tied 
closely to the existence of the Adler zero. Zheng et al. have reported similar 
results [159]. 

The E791 collaboration has also reported an observation of the k pole in L>+ — > 
K~7r~^7r~^ [160]. They parametrise the k width as F = Fop(s)/p(M^), where 
p(s) is Ktt phase space and find M = 797 ± 19 ± 43 MeV, Fq = 410 ± 43 ± 87 
MeV. The corresponding pole position is at M = 721—1292 MeV. The intensity 
calculated from their Breit-Wigner denominator is shown by the dashed curve 
on Fig. 73(c). Their k, is much narrower than from LASS data. 

BES II data for the k have been presented by Wu [161] in the channel J/^* — >■ 
K'^K~7r'^7r~. Strong K*{890) and 7^^0(1430) bands appear in the scatter plot of 
Fig. 74(a) and as peaks in (b). In (c) a peak due to p(770) is visible and in (d) 
peaks due to /^i(1270) and 7^'i(1410). The k appears as a broad low mass peak 
in the final state K*{890)K7i. It is shown in (f), where one K^n^ combination 
is chosen and then the mass distribution of the other K^tt^ combination is 
plotted; (e) shows the Dalitz plot of K*Ktt for this selection. 

In the analysis, it is necessary to show that the k signal does not arise from 
decays of A'i(1270) and -ft'i(1410), since all their decays give low mass Kti 
combinations. The broad k signal of Fig. 74(e) in fact arises partly from this 
source, but mostly from k, ir°(1430) and their interferences (and a small 
amount of ^^2(1430)). Two discrete solutions exist, both giving a k pole but 
with different widths. 

Fig. 75(a) shows that a fit without the k in the Ktt S-wave is bad. If the Adler 
zero is included in the formulae, there is a broad k pole with pole position 

M = (760 ± 20 ± 40) - i(420 ± 45 ± 60) MeV; the Ko{U30) pole position 
is (1433 ± 32) — i(181 ± 16) MeV. In this solution, there is rather strong 
destructive interference between k, and Kq^IASO). Contributions from k and 
Ko{lA30) alone are shown in Fig. 75(b) and further combinations in (c). The 
phase of the k is determinined accurately from interference between k and 
-^"0(1430). Results are shown as points on Fig. 75(d). They agree closely with 
the full curve, which shows the fit to LASS data. Note that in this solution, 
the coherent sum of k. and i^ro(1430) produces a low mass peak with a width 
of ~ 400 MeV. This narrower peak is similar to the solution found by E791. 
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Fig. 74. (a) The scatter plot of M{K^tt~) against M[K~tt~^). (b) The projection 
on to K^Ti^ mass; the histogram shows the fit. (c) The projection on to vrvr mass, 
(d) The projection on to K-ktt mass, (e) The Dahtz plot for events where one ivT^vr^ 
combination is in the mass range 892 it 100 MeV. (f) The mass projection of the 
second K^-k^ pair for the same selection as (e). 

The effect of the destructive interference is also to make the ii'o(1430) peak 
appear narrower than it really is. The width fitted to i^'o(1430) by E791 is 
low: 175 ± 12 ± 12 MeV. 

In the second type of solution, the Adler zero is omitted and a simple Breit- 
Wigner amplitude is fitted with F oc p{s). In this case, there is weak con- 



103 




Fig. 75. (a) The poor fit when the channel i^*(890)K is removed; (b) indi- 
vidual contributions from k (full histogram) and i^o(1430) (dashed); (c) the 
coherent sum of k + i^o(1430) (full histogram) and the coherent sum k + 
Ko(1430) + Ki(1270) + Ki(1410) (dashed); (d) the phase of the n amplitude (fuU 
curve) , compared with phases in individual slices of 100 MeV (points with errors). 

structive interference between k. and i^o(1430). The width fitted to the k, is 
smaller, 450-600 MeV, in agreement with the rather narrow k, peak in Fig. 
75(b). In order to fit this width of k, to LASS data it is necessary to add 
a strong background amplitude producing a negative phase; the sum of this 
repulsive contribution and the k generates the Adler zero [162]. 

Both approaches require a k pole. The difference between the two approaches 
reduces to one question. Should the Adler zero be included directly into the 
width of the k, or should it appear as a destructive interference between a 
conventional Breit-Wigner resonance and a background amplitude arising from 
the left-hand cut? Fitting both left and right-hand cuts together might resolve 
this question; but it is also possible that uncertainties about distant parts 
of the left-hand cut (how to include Regge poles) will obstruct this way of 
resolving the issue. Crystal Barrel data for pp — >■ KKn at rest find i^^o(1430) 
parameters M = 1415 ± 25 - i(165 ± 25) MeV. 
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11.2 Theoretical approaches 



Several theory groups have tackled the question of fitting Chiral Perturbation 
Theory (ChPT) to data. OUet, Oset, Pelaez and collaborators have published 
a series of papers where they use ChPT for the second order driving force and 
then unitarise the theory by including fourth order terms [163-168]. Some 
parameters of ChPT are not known accurately. They vary these and fit to 
data. A virtue of these calculations is that they also fit repulsive phase shifts 
in the I — 2 nn channel and the 1 — 3/2 Kn system. Poles are required for 
a, K, /o(980) and ao(980). A representative set of pole positions is quoted and 
fits to data are shown in a recent review by Pelaez [169]. Their a pole tends 
to come out around 440 — i215 MeV, slightly narrower and lower than values 
given above; their k is typically 750 — i226 MeV, again fairly narrow. 

Schechter and collaborators adopt an effective Lagrangian based on ChPT 
[170-173]. Their conclusion is that there is mixing between qq and a 4-quark 
system, similar to that proposed by Jaffe [174]. Beisert and Borasoy also adopt 
an effective Lagrangian approach [175] fitting experimental phase shifts but 
including r)' as an explicit degree of freedom. 

An approach based on meson exchanges was pioneered by the Jiilich group of 
Speth et al. [176-8]. Zou and collaborators adopt this approach fitting data 
with K- matrix unitarisation [179,180]. 

Morgan, Pennington and collaborators emphasise unitarity and analyticity as 
a general approach to fitting data [54,181,182]. Kaminski et al. hkewise fit 
data imposing constraints of analyticity and crossing symmetry [183]. The 
approach of Tornqvist is based on the linear a model, broken SU(3) and 
analyticity [184,185]. Fits by the Ishida group are reviewed in Ref. [186]. 

Important insight about the connection between mesons and their short-range 
content appears in the approach of Van Beveren and Rupp [187]. They write 
down two relativistic coupled Schodinger equations. One describes the qq sys- 
tem at short range. The second describes mesons at long range. There is a 
transition potential between them, based on the ^Pq model for transitions be- 
tween qq and 4 quarks. This models deconfinement of qq into mesons. Model 
parameters are fitted to observed spectra using a simple harmonic oscillator 
potential for confinement. Although the model is not precise, it gives con- 
siderable insight into the general picture. The Adlcr zero is built into the 
equations to a good approximation. The transition potential has the shape of 
the qq P-state and peaks at ~ 0.7 fm. 

The transition potential acts as a partially reffecting barrier, like the surface of 
a sheet of glass. The qq states sit in a confining potential, but can leak through 
the barrier. The interesting point is that a, k, /o(980) and ao(980) appear as 
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weakly bound states of the transition potential itself. They are mixed states of 
mesons and qq. Only S-wave mesons are bound by the transition potential. In 
this model, SU(3)-breaking appears directly through the masses of the mesons. 

I will risk offending everyone by summarising present problems. It is a mar- 
riage between theory and experiment: Lattice QCD calculations cannot yet 
contribute with any accuracy. Experimental approaches often fit only the 
physical region and ignore constraints from the left-hand cut; the work of 
Leutwyler et al. and Kaminski et al. are exceptions. All experimental fits have 
trouble with the opening of inelastic thresholds above 1 GeV. More work on 
understanding the 47r threshold could be profitable. 

Chiral Perturbation Theory is rigorously connected to QCD in principle, but 
in practice has a large number of parameters. The Adler zero clearly plays 
an essential role, but it is uncertain whether this should be built directly into 
resonance amplitudes or whether it appears as a cancellation between a con- 
ventional resonance and background from the left-hand cut. ChPT requires 
some unitarisation procedure and fitting parameters to data; its success verifies 
fundamental theoretical ideas. Meson exchange models suffer from uncertain- 
ties about heavy exchanges and Reggeisation. Analyticity and unitarity are 
sancrosanct, but in reality one is always battling with the effects of distant sin- 
gularities, not knowing quite what effects they have; Occam's razor is needed 
to get results. Theory does not know how to build in qq states and decays, 
except in the empirical approach of van Beveren and Rupp; that is where their 
scheme has much to offer, even if it is empirical. 

Despite rough edges between different approaches, the general consensus is in 
favour of a nonet of fight scalar mesons, with SU(3) breaking due to expficit 
mesons masses. The light mesons form a family of abnormal states resembling 
the deuteron and the virtual pp ^Sq state. 

11.3 /o(980) 

Recent BES II data on J/^P — 07r+7r~ and (pK^K^ provide a major improve- 
ment in the parameters of /o(980). It appears as a conspicous peak in Fig. 26, 
Section 5.5, interfering on its lower side with the a. The important point is 
that it is also observed clearly in decays to KK, Fig. 27. It is fitted with the 
Flatte formula: 

/ = 1/[M' - s - iiglM^^ + 9%KPis)KK)]- (38) 

The branching ratio between KK and tttt determines g^l Qi = 4.40 ± 0.25. In 
much earlier work, this ratio was unclear because of lack of information on 
KK. The value of g2 correlates rather strongly with M and gi because the 
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last term in eqn. (39) needs to be continued analytically below threshold so 

as to satisfy analyticity expUcitly. From BES data, M — 974 ± 4 MeV and 
gf = 142 ± 10 MeV. The paper gives the correlations. 

11.4 ao(980) 

Crystal Barrel data at rest provide accurate parameters for ao(980). Data on 
pp a;?77r° [188] determine the full-width at half-maximum very precisely: 
r = 54.12 ± 0.34 ± 0.12 MeV. This value is used as a constraint in fitting data 
on pp rj-K^Ti^. where the ao(980) is conspicuous [90]. On the Dahtz plot, 
the KK threshold appears as a very clear 'edge'. The intensity of the signal 
changes abruptly there. Prom the change in slope at this edge, the relative 
coupling to 777r and KK is determined rather accurately compared with most 
other experiments. The ao(980) in fitted with a Flatte formula like cqn. (39) 
with M = 999 ± 5 MeV, g^^ = 0.221 ± 0.020 GeV, g^K/g^^ = 1.16 ± 0.08. 

11.5 (p ^ 77r7r and ^yrjir 

Data for 0(1020) -frjn° and 77r°7r° from Novosibirsk [189-191] and the 
Kloe group at Daphne [192,193] are a topic of great interest. Achasov and 
Ivanchenko pointed out in 1989 that these data may throw light on ao(980) 
and /o(980) [194]. A vigorous debate has followed the publication of data; 
Boglione and Pennington [195] give an extensive list of references. 

1 have myself made unpublished fits to these data and wish to draw upon one 
important conclusion, which supports work of Achasov and collaborators. The 
fits are made using parameters of /o(980) and ao(980) given in the previous 
subsections; formulae are taken from the KK-\oop model of Achasov and 
equations of Achasov and Gubin [196]. Apart from the absolute normalisation, 
there are no free parameters. The results of the fit to Kloe data for the rin'-f 
final states arc shown in Fig. 76(a) and (b). The good fit immediately verifies 
the KK loop as the correct model for the process. This loop gives rise to a 
violent energy dependent factor |/(a,6)p shown in (d). This multiphes the 
factor for the Ml transition; the product of all energy dependent factors in 
shown in (c). This leaves httle doubt that the basic mechanism is understood. 

The situation is not quite so good for (f) — > 77r°7r°. Pig. 77 shows the fit. 

However, there are two complications which affect the mass range below 0.6 
GeV (not important for present conclusion). The first is that the final state 
pn'^, p 7r°7 contributes there and is poorly known in normalisation. The 
second is that there may be a small a contribution there. The fit to the /o(980) 
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Fig. 76. Kloe data for 4> 'yrjir^ for (a) ry 77, (b) rj 7r+7r~7r'^; full curves 
show the fit. The dashed curve in (a) illustrates uncertainties from a possible form 
factor and the Adler zero; (c) the energy dependent factor appearing in dV/dm, but 
omitting |/(o,6)p, (d) \I{a,b)\^ itself. 
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Fig. 77. Kloe data for (j) 77r tt ; the full curve shows the fit and the dashed curve 
the /o(980) intensity alone. 

peak using BES parameters looks slightly narrow, but can be optimised by 
increasing g'^^ slightly to 160 MeV. With those mild reservations, the fit is 
again good to the /o(980) peak and verifies the i^^-loop model. 

Achasov argues that these data reveal the origin of ao(980) and /o(980) [197]. 
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Fig. 78. Momentum-space wave functions of /o(980), full curve and ao(980), dashed. 

The transition (p ao(980) is OZI- violating, since ss must change into nn. 
Physically, the argument is that ss states like (j) and /2(1525) do not decay 
readily to non-strange mesons. His conclusion is that ao(980) must be a state 
such as Jaffe [174] has proposed, made up of (ussu — dssd)/\/2. The u and d 
quarks in the i^^-loop appear in the /o(980) and ao(980) final states. 

Jaffe's state involves 4 quarks with a small spatial extent; the binding en- 
ergy arises from pairing of quarks. Pairs of q and q form an SU(3) flavour 3 
and 3 configurations; these combine to a fiavour singlet. Note that all these 
configurations are confined to rather short ranges. 

A KK molecule contains the same quarks, but configured into KK. It is 
not obvious how to separate these two possibilities. The /o(980) and ao(980) 
have small binding energies. They must therefore be surrounded by long-range 
clouds of kaons. For a given binding energy B, the wave function of the kaon 
pair is just oc (1/r) exp(— //r), where fi oc y/MB; M is the reduced mass. 
Figs. 78 and 79 shows results obtained by weighting \t,b\'^ with the line-shape 
of /o(980) or ao(980). Wave functions of ao(980) and /o(980) are remarkably 
similar. This is illustrated in momentum space k in Fig. 78. It need not be so. 
The Flatte formulae describing them have very different parameters at face 
value. Nonetheless, the wave functions of /o(980) and ao(980) are within 10 
percent of one another for all k. That suggests they have a common character. 

Tornqvist [18] points out that both ao(980) and /o(980) must contain a KK 
component. His eqn. (15) includes mesonic channels A^Bf. 

^ \qq>+Zr[-{d/ds)Re H,(g)] > 

l-Eiid/ds)ReUi{s) ' ^ ' 

Here Re H(s) = glcK\J ^''^'kI ^ ~ 1 ^r s < Am%;. This formula is easily eval- 
uated to find the KK components in ao(980) and /o(980) as functions of s. 
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At the KK threshold, these components are of course 100%. But they are 
> 40% for all masses above 900 MeV. This is illustrated in Fig. 79 by the 
dotted curve. This figure also shows as the full curve the line-shape of /o(980) 
V. mass. The RMS radius of the kaonic wave function is shown by the dashed 
curve in Fig. 79, scaled down by a factor 10 for display purposes. When the 
intensity of /o(980) is above half-height, the RMS radius is > 1.7 fm. 

The photon in 4> radiative decays couples to kaons dominantly at large r, 
where the El matrix element has a large dipole moment e\r\. Note too that 
the kaons couple with charge 1, unlike the quarks which have charges 1/3 or 
—2/3. Achasov argues that KK 'molecules' cannot explain the observed large 
rates of (f) radiative dedcays. However, a quantitative assessment of the rates 
is needed, taking into account the large radii of the kaonic clouds. It is hard 
to believe that the matrix element for radiative decays can be dominated 
by small radii. 

However, the matrix element for KK /o(980) or ao(980) will be concen- 
trated at small r. This matrix element may go via a short-range interaction at 
the quark level. Or it may go through ^-channel mesons exchange; these will 
be dominated by exchange of K*[890). This exchange will be responsible for 
TTTT KK. It is then 'natural' that ao(980) and /o(980) lie close to the KK 
threshold. In Ref. [150], the matrix element for nn tttt was found from the 
Fourier transform of the matrix element in momentum space. Although this 
matrix element is uncertain above 1 GeV, as shown by the dashed curve in 
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Fig. 79. The line-shape of /o(980), full curve; the KK fraction in the wave function, 
dotted curve; and 0.1*RMS radius, dashed curve. 
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Fig. 80. Matrix elements F for tttt — tttt (a) in momentum space, (b) F(r) 

and (c) r^F^{r). The full curves show the optimum fit to elastic data and dashed 
curves show results with a stronger exponential cut-off in s above 1 GeV^. 

Fig. 80(a), F{r) is definitely concentrated at a radius < 0.4 fm. It appears that 
binding of /o(980) and ao(980) could come cither from interactions between 
quarks or from meson exchanges. These states must have large inert kaonic 
clouds, as Tornqvist emphasises. 

Kalashnikova and collaborators [198] re-examine an old idea of Weinberg to 
decide whether the ao(980) and /o(980) are composites of elementary particles 
or are 'bare' states. The choice between these possibilities rests upon whether 
the resonances have a single pole close to the KK threshold or a pair of 
poles on different sheets. A single pole favours a composite, like the deuteron. 
Morgan made this point in a different way in 1992 [199]. The BES data give 
a clear answer for /o(980). There is a single nearby second-sheet pole at M = 
0.999 -i0.014 GeV and a distant third-sheet pole at M = 0.828 -iO.345 GeV. 
For ao(980), the situation is not so clear-cut. There is a second-sheet pole at 
1.032 - i0.085 GeV and a third-sheet pole at M = 0.982 - i0.204 GeV. 
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11.6 L>,(2317) 



The discovery of the narrow Z^s(2317) by Babar [200] and Cleo [201] has gener- 
ated a flurry of interest. The predicted mass for a O"*" cs state in typical quark 
potential models is 2.48-2.49 GeV, above the DK threshold at 2360 MeV. 
Van Beveren and Rupp [202] argue that the Ds(2317) is a cousin of the light 
scalar nonet associated with the DK threshold. They present a calculation 
following an earlier one of ao(980) [203]. 

There, they model the attraction in the 1=1 KK S-wave, ignoring coupling 
to rj-K. Fig. 81 shows the movement of the pole as the coupling constant to 
KK is varied. A virtual state appears first at 0.76 GeV. As the coupling to 
KK increases, the pole migrates along the lower side of the real axis to the 
KK threshold. For still larger coupling, it becomes a bound state and slowly 
moves away from the KK threshold on the upper side of the cut. The pole 
remains in the vicinity of the KK threshold for quite a wide range of coupling 
constants. This leaves little doubt that the KK threshold plays a strong role 
in the dynamics. If the mass of the kaon were to change, the pole would move 
with it. This does not necessarily mean that ao(980) is a bound-state of KK. 
The attraction in this channel may arise from a variety of processes generating 
attraction, either at the meson level or at the quark level. 

Likewise, the Ds(2317) appears as either a virtual or bound state of DK with 
some cs component mixed in. The regular cs state is predicted to move up in 
mass above the mass of simple quark models through level repulsion. A similar 
phenomena near the BK threshold is of course predicted. However, there are 
many other competing ideas to explain Ds(2317) and other narrow states. 

The (7 and k appears close to tttt and K-k thresholds, but are forced away 
I^e(E) 0.70 0.80 0.90 1.00 GeV 



-1.0 GeV 



-2.0 GeV 



Im(E) 



Fig. 81. Prom Ref. [203]. The movement of the ao(980) pole is shown as a function 
of the coupling constant for the 1=1 KK S-wave. Values of coupling constants are 
indicated. The inset shows a magnified view near the KK threshold. 
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from those thresholds by the Adler zeros. It is good for us that this happens. 

If they were bound states below threshold, the Universe as we know it would 
be totally different because of the long-range o exchange binding nuclei. Van 
Beveren and Rupp predict [202] that there will be a Dtt 0^ state not far above 
the D-K threshold. This could be the ^0(2290) reported by Belle with a width 
of 305 MeV [204]. 



12 Future Prospects 

It is possible to check experimentally whether /o(1710) is a ground-state qq 
state or a radial excitation. This may be done by checking its production rate 
in /o(1710)7r. Van Beveren, Rupp and Scadron argue [205] that /o(980) 

is dominantly ss because of its absolute production rate. This is calculated 
assuming that the decays by c — > followed by turning into 

Tr"*". Their argument could be perturbed by other diagrams, but these will 
only change the rate by small factors. Implicit in their calculation is that 
there is good overlap between the radial wave functions of Ds and the radial 
wave function of the /q. This is true for the mesonic ground-state. But if the 
/o(1710) is a radially excited ss state, as Anisovich and Sarantsev would have 
us beheve, the decay rate will be suppressed by at least a factor 100. 

E852 and VES have data which could locate the missing qq states in the mass 
range 1600-1800 MeV. The 7 = hybrid partner of 7ri(1600) is also missing. 

CLEO C is already taking data on data on J/^ and \E'' will follow during 
the next 5 years with better angular acceptance than BES II and excellent 
7 detection. A high priority should be looking for the 2''" glueball in ■0" — > 
77(477) as suggested above. Decays J/^ — > ^r\r\ and 700 are important for 
full information on O"*", 2"*" and 0~ glueballs and sorting out mixing schemces 
with qq. A search for /o(2100) KK and 7777 would be simple and valuable. 
Radiative decays of are almost unknown at present; the node in its radial 
wave function will throw light on wave functions of resonances produced in 
^' decays. High quahty data on x states should also appear. This looks an 
exciting program of physics over the next few years. 

The way forward with 1 qq states is to study diffractive dissociation of the 
photon. 

What would complete the picture of qq states and make it more secure is to 
rerun the Crystal Barrel experiment using a frozen spin polarised target. Study 
of 7r°77 and 37r° would remove ambiguities and tighten the analysis for 7 = 1, 
C — Study of a;7r° and a;77 channels would likewise improve and hopefully 
complete 7 = 1 and 0, C = — 1 states. Such an experiment is technically 
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straightforward. Data are needed right down to 360 MeV/c to cover the tower 

of resonances around 2 GeV. This is a possible program using the new p ring 
proposed for GSI. It could bring the long program of work on light non-strange 
mesons to a definitive conclusion. 
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